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1.  Introduction. 
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The  gonl  of  this  project  has  been  to  study  the  estimation  and  control  of  elastic  systems  com¬ 
posed  of  beams  and  plates.  This  general  goal  has  served  as  a  guide  for  our  research  over  the  last 
several  years.  We  have  made  substantial  progress  in  our  investigation  of  topics  within  this  and 
related  areas. 

For  control  we  have  studied  two  basic  problems.  The  first  is  that  of  locating  the  optimal 
placement  of  controllers  on  a  beam  or  a  plate  in  the  case  of  either  static  or  dynamic  models. 
The  second  involves  the  control  of  general  three  space  dimensional  elastic  models  that  incorporate 
nonlinear  friction  and  contact  laws  in  their  boundary  conditions. 

In  estimation  and  identification  we  have  obtained  results  concerning  the  estimation  of  elastic 
coefficients  in  beams  and  plates  for  static  and  dynamic  models.  Further,  we  have  investigated 
the  estimation  of  certain  damping  terms  arising  in  plate  models.  Finally,  in  related  work  we 
have  considered  the  determination  of  diffusion  parameters  in  parabolic  equations,  and  conditions 
for  bgectivity  of  the  parameter  to  state  mapping  for  discrete  approximations  of  certain  elliptic 
boundary  value  problems. 

Our  work  has  sought  generally  to  explore  the  theoretical  analysis  of  the  distributed  models, 
algorithms  for  numerical  implementation,  and  an  approximation  theory  relating  the  two. 


2.  Control  Problems 

Let  n  be  an  open  bounded  domain  in  K’  with  a  Lipschitz  boundary  F  which  is  occupied  by  a 
plate.  We  consider  dynamic  models  of  the  form 


(2.1) 

where 


and 


utt  flut  -i-  Au  =  in  n  X  (0,T) 


Bifi  =  A(6oA^)  -  V '  {biV<p)  +  bfip 


A(p  =  A(aoAy7). 
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We  have  assumed  simply  supported  boundary  conditions  although  our  results  carry  over  to  other 
boundary  conditions  as  well.  Coefficients  ao,h,^i,^3  L°°(n)  satisfy  assumptions  such  as  uq  > 

u  >  0,  bi  >  0  while  #  is  taken  to  belong  to  Defining  the  control  functional  J  to  be  for 

example, 

(2.2)  J{0)=  f  f  {u{x,t\0)-z{x,t)Ydxdt+  f  {eJl  +  eofi^)dt 

Jo  Jn  Jo 

with  si.ea  >  0  and  Ci  +  cj  >  0,  we  consider  the  problem:  find  in  H^{0,T){L^{0,T))  such  that 

(2.3)  Ji0o)  =  mfimum{J(/3)  :  ^  €  ff‘(0,r)(L»(0,r))} 

The  actuator  is  modeled  by  means  of  the  term  In  fact,  $  may  be  chosen  as  a  Dirac  measure 
6x„  fo  xq  €  n  or  as  a  convolution  with  a  function  ^  =  ^(*  -  xo)  for  lo  G  D  with  ^  G  ^*(n)  (or 
smoother)  and  aupp^  C  {i :  |i  -  zo|  <  7}  where  'j  <  dist{xo,T)- 

The  existence  of  a  unique  solution  to  problem  (2.1)-(2.3)  follows  from  standard  arguments  [2]  . 
Of  interest  in  this  investigation  is  the  dependence  of  the  optimal  control  ^0  and  the  functional 
J{0o)  upon  In  particular,  for  $  =  ^*0  or  $  =  ^(’  -  iq)  the  uniqueness  of  the  optimal  control 
implies  that  the  mappings 

Xq  •-*  ;3(xo)  and  xq  •-»  y^xo)  =  J{0o{xo)) 

are  well  defined.  One  may  now  consider  properties  of  these  mappings.  Indeed,  if  continuity  can 
be  established,  then  given  any  closed  subset  F  of  H  there  is  a  best  point  x  of  control  among  those 
points  of  F  in  the  sense  j(*o)  <  {j(*)  :  *  €  F}.  If  differentiability  can  be  established,  then  it  may 
be  used  to  construct  algorithms  to  find  the  optimal  control  location.  Roughly,  we  show  that  if  the 
difference  quotients 

=  («(xo  +  €h)  -  ♦(*o))/e 

converge  in  to  an  element  4',  then  the  mappings 

lo  *-*  0{xo)  and  Xo  y(*o) 

are  differentiable,  cf.  [3].  Furthermore,  the  derivative  of  ^(xq)  is  shown  itself  to  be  a  solution 
of  an  optimal  control  problem.  In  fact  the  derivatives  of  0{xo)  may  be  associated  with  a  sequence 
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of  optimal  control  problema.  The  length  of  the  sequence  is  dependent  upon  the  differentiability 
of  These  properties  are  used  to  compute  optimal  points  of  control  in  several  test  cases  using 
descent  algorithms  for  beams  and  plates  in  [9,10].  Similar  consideations  are  given  to  the  static  case 
in  [8]. 

During  the  grant  period  we  also  studied  the  control  of  certain  models  of  three  dimensional 
elastic  bodies.  These  models  are  in  the  form  of  hyperbolic  variational  inequalities  and  incorporate 
nonlinear  normal  and  frictional  interface  laws  in  through  their  boundary  conditions.  They  were 
developed  by  J.T.  Oden  and  his  coworkers  [5,6]  and  seem  to  produce  results  that  agree  well  with 
experimental  observations  [5].  Since  these  models  are  rather  complicated,  we  refer  the  reader  to 
Appendix  A  for  their  formulation. 

We  have  obtained  results  on  optimal  control  by  means  of  distributed  body  forces  and  boundary 
forces.  The  existence  of  optimal  controls  we  have  proved  for  the  variational  inequalities  directly. 
Since  certain  functionals  within  the  variational  inequality  are  not  differentiable,  we  consider  regu¬ 
larizations  of  the  functionals  and  of  Sobolev  type.  For  the  resulting  problems  we  obtain  optimality 
and  regularity  conditions  for  solutions.  Using  these  results,  we  proceed  to  obtain  convergence 
results  for  finite  element  approximatons  to  the  regularized  problems  and  iterated  limit  theorems 
relating  convergence  to  the  optimal  controls  of  the  variational  inequedity  see  Appendix  A. 


3.  Identification  and  Elstimation. 

Let  n  be  an  open  bounded  domain  in  S’  with  a  Liptschitz  boundary  F  and  let 


(3.1)(a)  Au  =  A(aoAu)  —  V  •  (oiVu)  +  aju 

or 


(3.1){6)  Au  =  <rA(ooAu)  +  (1  -  <r)((ao«.»)s»  +  +  2(oou,,),,) 

and 


(3.2) 


tttt  +  (A(6oA)  —  V  •  (6iV)  +  h^)ut  +  At*  —  /  in  0  x  (0,  T) 


accompanied  with  appropriate  boundary  and  initial  conditions.  Equations  (3.1)-(3.2)  along  with 
appropriate  initial  and  boundary  conditions  prescribe  a  mathematical  model  of  a  thin  plate  com¬ 
monly  called  the  Kirchhoff  plate  model  [7].  Also,  terms  are  included  in  (3.2)  to  model  damping. 
The  choices  of  these  terms  are  motivated  from  analogous  beam  models  [4]. 

The  estimation  problem  may  be  formulated  as  follows.  Given  information  z,  which  we  view 
as  a  member  of  an  observation  space  Z  (Hilbert  space),  we  seek  to  determine  one  or  more  of  the 
parameters  q  =  {<^,6^  :  t  =  0, 1,2}  from  within  a  prescribed  subset  Qad  of  a  Hilbert  space  Q  that 
produces  a  solution  u{q)  of  (3.2)  that  “matches”  z.  One  approach  to  this  problem  sometimes  called 
the  (regularised)  output  least  square  method  is  formulated  as  a  minimization  problem: 

Find  90  €  Qad  such  that 


(3.3)(a)  J{qo)  =  infimum{J{q)  :  q  G  Qad} 

where  for  ^  >  0 


(3.3)(6)  J{<l)  =  \\Cuiq)-z\\l+mh 

and  C  is  the  observation  operator  taking  a  solution  0(9)  to  its  image  Cu{q)  in  the  observation 
space  Z.  The  set  Qad  is  specified  by  constraints  that  assure  problem  (3.2)  is  well-posed  and  that 
guarantee  sufficient  compactness  to  imply  existence  of  a  solution  to  (3.3). 

To  solve  these  problems  numerically  one  must  approximate  the  boundary  value  problems  by 
a  system  of  equations  and  approximate  parameters  in  some  consistent  manner.  The  adaptation  of 
the  constraints  to  the  discrete  problem  is  then  necessary.  We  have  investigated  this  problem  for 
static  beams  and  certain  static  plates  in  [11,12]  and  [13].  These  papers  consider  penalization  and 
regularization  techniques  to  incorporate  constraints.  Useful  here  is  an  analysis  of  the  regularity  of 
the  optimal  estimators  of  these  problems  and  the  use  of  linear  splines  to  approximate  parameters. 
Since  linear  splinos  preserve  pointwise  bounds,  such  pointwise  L°°-constraints  easily  carry  over  to 
finite  dimensional  problems.  Use  of  this  technique  however  restricts  the  class  of  problems  that  may 
be  considered  for  plates  to  those  in  which  the  domain  O  is  a  rectangle  and  the  parameters  are 
expressible  as  a  tensor  product  of  spaces.  This  second  restriction  is  not  so  great  since  such 
spaces  are  dense  in  C*’(n).  However,  the  theory  is  still  essentially  a  one  dimensional  theory. 

For  the  general  case  we  have  studied  a  two  dimensional  theory  [13,14]  in  which  Qad  >•  assumed 
to  be  in  l?’(n).  For  this  case  there  typically  is  not  a  finite  dimensional  ^>proximating  basis 
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in  H*{0)  that  preserves  pointwise  bounds  and  remains  in  H*{n).  The  regularity  is  determined 
by  appropriately  applying  a  generalized  version  of  the  Kuhn-Tucker  Theorem  and  analyzing  the 
regularity  of  the  solutions  of  the  resulting  vai.ational  equations.  Using  these  properties  we  may 
appropriately  define  finite  dimensional  sets  of  parameters  to  obtain  convergence  theorems. 

For  example  in  [14]  we  consider  the  problem 


utt  +  Au  =  f  in  n  X  (0,  T) 
u=^=0  on  Fx  (0,r) 

(3.4) 

u(0)  =  uo  in  n 
ut(0)  =  ui  in  n 

where  A  satisfies  (3.1)  with  ai  =  oj  =  0  and  a  is  Poisson’s  constant  in  (0,  j).  Of  great  importance 
in  this  analysis  is  the  determination  of  regularity  properties  of  the  optimal  estimator.  If  F  is 
and 

/  €  L*(0,r))  n  ff*(0,T;  ^’*(0)) 

Uo  €  D(A) 

«x  €  ffo’(n) 

and 


z€L»(0,r;L>(n)), 


consider  problem  (3.3)  with 


=  {a  €  ’(n) :  a(x)  >  »^  >  0  in  fi}. 

Of  course,  other  functionals  J  may  be  considered  as  weU  by  matching  at  time  T  or  at  points  (  within 
n.  We  obtain  regularity  results  indicating  that  the  optimal  parameter  oq  belongs  to  for 

3  €  [0,1).  This  result  is  then  useful  in  formulating  finite  dimensional  approximating  problems 
based  on  Galerkin  approximation  of  (3.4)  and  in  specifying  finite  element  subspaces  in  H*{Q)  to 
approximate  the  coefidents  a. 

Numerical  examples  have  been  studied  for  this  problem  in  Q  =  (0, 1)  x  (0, 1)  and  T  =  0.5. 
Tensor  products  of  cubic  B-splines  on  a  regular  mesh  adjusted  to  boundary  conditions  are  used  to 
approximate  u.  Again,  tensor  products  of  cubic  B^splines  on  regular  mesh  are  used  to  approximate 


'v*  ',1 


a.  The  reeulting  system  of  ordinary  differential  equations  is  approximated  using  a  symmetric  time 
differencing  technique.  The  following  examples  use  49  basis  elements  to  approximate  the  state  u 
and  25  to  approximate  a.  The  parameter  0  equals  IQ~*,  and  the  test  coefficient  is 


CLut  =  1  +  (i*  +  y’)/2 


with  observation 


2(i.y,0  =  16i*(l  -  x)*y*{l  -  y)*  cos(t) 

The  forcing  function  /  is  generated  using  ot,t  and  z.  Using  a  quasi-Newton  method  for  minimization 
with  initial  guess  a  =  .25,  an  approximation  for  at,t  with  relative  error  of  2.7%  is  obtained  after 


7  iterations. 

For  a  similar  case  with 


=  l/(4((i  -  i)*  +  (y  -  i)*)  +  1) 

an  ^proximation  for  <it*t  with  relative  error  of  2.5%  is  obtained  after  5  iterations. 

We  have  analyzed  similar  problems  for  the  estimation  of  damping  terms  as  in  (3.2)  with 
6i  =  6,  =  0  and  have  obtained  regularity  and  approximation  results  [15].  With  a  test  coefficient 

»,«  =  l/(4((x-l)’  +  (y-l)')  +  l) 

with  initial  guess  having  relative  error  16%  after  4  iterations  we  obtain  a  relative  L*  error  of 
3.8%  with  10-“. 

In  other  work  we  have  studied  the  bijectivity  or  identifiability  of  the  parameter-to-state  map¬ 
pings  for  finite  dimensional  version  of  elliptic  boundary  value  problems  [ij.  In  other  work  with 
K.  Kunisch  we  have  studied  estimation  of  time  and  spacially  dependent  diffusion  coefficients  in 
parabolic  equations  see  Appendix  B.  In  this  paper  we  obtain  regularity  properties  of  the  estimated 
diffusion  coefficient  and  also  determine  sufficient  conditions  for  strict  complementarity. 
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Abstract.  In  this  paper,  we  study  the  control  of  the  motion  of  viscoelastic  bodies  resisted  by 
contact  and  frictional  forces.  The  contact  and  friction  effects  are  modeled  by  new  interface  laws 
that  characterize  friction  phenomena  on  dry  rough  surfaces.  The  system  under  study  is  an  arbitrary 
body  composed  of  a  linearly  viscoelastic  material  which  is  subjected  to  body  forces  and  tractions 
and  which  is  in  contact  with  a  rough  surface  of  a  neighboring  body.  Under  mild  conditions  on  the 
data,  we  are  able  to  prove  the  existence  of  an  optimal  control  of  motion  of  such  systems  in  L- 
(O.T;V)  through  a  program  of  applied  body  and  surface  forces. 
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1.  INTRODUCTION 

The  dynamics  and  control  of  distributed  systems  subjected  to  hricdonal  resistive  forces  has 
been  an  open  problem  in  partial  differential  equations  for  many  years,  largely  owing  to  the  absence 
of  a  meaningful  existence  theory  for  dynamical  systems  with  friction.  The  mathematical  difficulties 
associated  with  Coulomb  Mction  in  problems  of  elastodynamics  were  pointed  out  in  1972  by 
Duvaut  and  Lions  (see  [1])  and  has  been  the  subject  of  much  study  in  the  intervening  years.  The 
finite-dimensional  case  has  proved  to  be  somewhat  more  tractable,  and  conditions  fc^  the  existence 
of  solutions  of  discrete  dynamical  problems  with  friction  have  been  recently  reported  by  Lotstedt 
[2]  and  Jean  and  Pratt  [3].  The  complexity  of  the  problem  can  be  appreciated  by  reviewing  the 
work  on  dynamics  of  systems  with  frictionless  contracts  by,  for  example,  Schatzman  [4],  Carrero 
and  Pascal!  [5],  Ldtstedt  [6],  Amerio  and  Prouse  (7],  Amerio  [81,  and  others. 

In  recent  papers,  Oden  and  Martins  [9,10]  pointed  out  that  one  of  the  principal  sources  of 
mathematical  difficulty  was  the  definition  of  frictional  stresses  on  the  contact  surfaces  characterized 
by  Coulomb's  law.  However,  an  overwhelming  volume  of  experimental  data  accumulated  over  a 
half  century  suggests  that  this  law  is  inadequate  for  modeling  actual  contacts  and  resistive  forces  on 
deformable  bodies.  By  characterizing  the  actual  normal  compliance  of  elastic  interfaces,  a 
constitutive  equation  for  an  interface  can  be  developed  which  yields  results  in  agreement  with  a 
sizable  collection  of  experimental  results  on  static  and  dynamic  friction  [9].  Moreover,  the  use  of 
such  interface  constitutive  laws  in  mathematical  models  of  elastostatics,  elastodynamics,  and 
viscoelastodynamics  problems  with  friction  produces  a  tractable  theory:  results  on  the  existence 
and  local  uniqueness  of  solutions  to  static  problems  in  elasticity  with  the  Oden-Martins  [9] 
nonlinear  friction  laws  were  recently  established  by  Rabier  et  al.[  1 1  ]  and  of  dynamic  problems  on 
elasticity  and  viscoelasticity  by  Martins  and  Oden  (10).  These  new  theories  and  results  set  the  stage 
for  the  study  of  the  optimal  control  of  such  systems,  taken  up  in  the  present  paper. 

In  the  present  study,  we  establish  the  existence  and  uniqueness  of  a  class  of  optimal  controls 
for  a  broad  class  of  problems  in  the  dynamics  of  linear  viscoelastic  bodies  with  contact  and  friction 
laws  of  the  type  introduced  by  Oden  and  Martins  [9].  We  are  able  to  show,  under  very  mild 

restrictions  on  the  data,  that  for  a  control  of  the  form  ^(t)H,  H  being  a  functional  characterizing  all 

external  forces  on  the  system,  a  |3e  H*  (0,T)  can  be  found  which  minimizes  a  functional  relating 
the  L^-distance  between  solutions  to  the  viscoelastodynamics  problem  and  an  arbitrary  target 


motion  z(*,t)  over  a  time  interval  [O.T]. 


Following  this  Introduction,  we  introduce  in  Section  2  notations  and  record  preliminary 
results  on  the  formulation  of  the  contact  problem  in  viscoelastodynamics.  In  Section  3,  a  weak  or 
variadonal  formulation  of  the  problem  is  presented  and  an  existence  result  from  [10]  is  recorded. 
The  optimal  control  problem  is  introduced  in  Section  4  and  the  major  results  are  established  in 
Section  5. 

We  believe  that  these  results  are  the  first  of  this  type  on  optimal  control  of 
distributed-parameter  systems  in  contact  with  friction. 

2.  PRELIMINARIES 

We  begin  by  considering  a  metallic  body,  the  interior  of  which  is  a  bounded  open  domain  Q. 
in  IR^  (N  =  2  or  3).  The  boundary  F  of  Q  is  smooth  (e.g.,  at  least  Lipschitz  continuous)  contains 
open  subsets  Fq,  Fp,  and  F^^.  such  that  F  =  Fp  u  Fp  u  F^,,  meas  =  0,  a  6  {D,F,C}. 

Material  particles  (points)  in  Q  with  canesian  coordinates  Xj;  .  I  :S  i  ^  N,  are  denoted  x,  the 

volume  measure  by  dx,  and  the  surface  area  measure  by  ds,  and  a  unit  exterior  vector  normal  to  F 
by  n. 

We  shall  assume  that  the  hxxiy  is  composed  of  a  linear  viscoelastic  material,  the  mechanical 
response  of  which  is  characterized  by  the  constitutive  law, 

~  ^ijkl  ^k.l  ^ijki  ^k.l 

1  ^  i,j,k.l  <  N  (2.1) 

where  (»  Ojj)  are  the  canesian  components  of  the  Cauchy  stress  tensor,  Ejjy  =  Ejjy  (x)  and 

Cijy  =  Cjju  (x)  are  the  arrays  of  elastic  and  viscoelastic  parameters  at  point  x  e  Q.  U|j  =  Uj^  (x,t) 

are  the  components  of  displacement  at  x  at  time  t6[0,T]  IR,  U|j  |  =  9uj^dx|  and  (*)  denotes 
differentiation  with  respect  to  time.  The  body  is  subjected  to  body  forces  of  intensity  b  per  unit 

volume,  to  surface  tractions  t  on  Fp,  and  is  possibly  in  contact  with  a  moving  neighboring  body 


on  the  contact  surface  F^,  which  moves  relative  u>  Q  at  a  velocity  U-p  tangent  to  F^,.  The  partial 

differential  equations,  inequalities,  and  condidons  governing  the  behavior  of  the  body  are  listed  as 
follows: 

Linear  Momentum 

cr.j  (u  ,  ij).j  +  bj  =  puj  in  Q  X  (0,T)  (2.2a) 

where  (u  ,  u)  is  the  expression  on  the  right-hand  side  of  (2.1). 


Rniin/lnrv  C.nnditinnx 


Uj  =  0  on  Fpj  X  (0,T) 
Ojj(u,u)nj  =  tj  on  Fp  x  (0,T) 


Contact  Interface  Conditions 


On  (u.u)  =  -  C„h  (u/  " 
o-p(u.u)=  0 

lo-p  (u.u)l  <  Cph  (u^r  ^ 
and 

IOt  (u,u)l  <  Cjh  (u^f  ^ 


“n>8 


Wj  =  0 


lap(u,u)l  =  Cjh(Un) 


niT 


=>  there  exists  X  >  0 
such  that 


Wj  =  -  X  o-p  (u.u) 
on  F(-  X  (O.T) 


(2.2b) 

(2.2c) 


(2.2d) 


(2.2e) 


r^itin!  r^rnHitums 


U  (X,0)  -  Uq  (x) 
ti  (x,0)  =  Uj  (x) 


x£  n 


(2.20 


In  (2.2a) ,  p  »  p(x)  is  the  mass  density  and  in  (2.2d)  and  (2.2e), 

(u,u)  a  the  nonnal  stress  (contact  stress)  on 

a  Cy(u.u)ninj 

h  (u„)  a  the  approach  of  the  two  contact  surfaces 
*  (Un  -  g)+ 

with  Ufl  a  ujnj  the  nonnal  component  of  displacement,  g  the  initial  gap  between  surfaces,  and 
(*)+  the  positive  part  of  the  indicated  argument  (4>^  =  max  (0,0)). 

CT.J.J  (u.ii)  a  the  tangential  (frictional)  stress  component  on  , 

a  0-- (u,u)  nj  -  njC^  (u.u) 


Wj  a  the  relative  (slip)  velocity  of  the  contact  surfaces. 


lly  -  U*^ 


where  u-j-j  =  uj  -  nj  u„. 


I 


In  (2.2d)  and  (2.2e).  Cg.  nig.  Cj->  ^  material  constants  characterizing  the  mechanical 

response  of  the  interface.  Equation  (2.  d)  is  the  power-law  constitutive  equation  for  the  nortnal 
compliance  of  the  interface;  (2.2e)  defines  the  corresponding  tangendal  response  characteristics. 
Justification  for  the  use  of  such  interface  laws  is  given  in  ODEN  and  MARTINS  [9].  These  laws 
correspond  to  a  generalized  Coulomb  law  in  which  the  coefficient  of  friction  is  dependent  on  the 

deformation  and  is  of  the  fomi,  |j.=  Clag(u.u)l®  with  a  a(n>j/mg)-  1  and  C  =  Cj/Cg"''^"'T 


3.  A  VARUTIONAL  FORMULA 


We  now  record  a  weak  formulation  of  Problem  (2.2),  in  which  regularity  requirments  on  u 
and  on  the  data  are  relaxed.  We  introduce  the  spaces  of  admissible  functions. 

V  =  (v  =  (vj.  V2, . . ..  vj^)  €  (H*(Q)*^I  v  =  0  a.e.  onTj^j 

(3.1a) 

H  =  (L2(a))N 

(3.1b) 

V'  =  topological  dual  of  V 

(3.1c) 

The  space  V  is  equipped  with  the  usual  norm. 

"v  II  =  {  In  (VjjVj.  +  vjvj)  dx}  1/2 

and  (*,‘)  and  1  ’  1  denote  the  L^-inner  product  and  norm  on  H,  and  <*,'>  denotes  duality  pairing  on 
V'  X  V.  Throughout,  we  assume  for  simplicity  that 

p(x)  *  1  and  meas  (Tq)  >  0. 

(3.2) 

In  (3.1a),  V  on  the  boundary  is  interpreted  in  the  usual  sense  of  traces  of  H*(Q)  -  functions. 


The  coefficients  in  (2.1)  are  assunxd  to  satisfy  the  following  conditions: 

^iikl*  Qikl  ®  L^’CQ) 


Hjld  ^ 

^ijkl  “  ^jilk  ~  ^ijlk  ~  ^klij 


^ijkl  “  ^jikl  =  Qjik  =  ^klij 


a.e.  X€  (2 


3  o^,  >  0  such  that 

for  every  symmetric  matrix  Ajj  e 


Whenever  (3.3)  hold,  the  bilinear  forms  a:  V  x  V  IR  and  c:  V  x  V  IR  (representing  the 
internal  virtual  work  of  the  stress  o--)  defined  by 


a(v.w)  =  Jjj  E.jyv^lWjjdx 


=  In  Cjjyv^ 


(3.4a) 


(3.4b) 


for  v,w  €  V,  are  continuous  and  V-clliptic,  i.c,  there  exist  positive  constants  M£,  M^,  Og.ct^;; 
such  that  for  any  w,  v  in  V, 


la(w,v)l  i  Mg  llwll  llvll  ,  lc(w,v)l  $  M^;.  Ilwllvll 


a(v,v)  k  Og  llvll^  .  c(v,v)  >  llvll^ 


(3.5a) 


(3.5b) 


s 

A 

4 

•  I 

5 


>1 

•i 

I 

f 


We  also  assome  diat 


1  ^ 


<+»•  if  N  =  2 
^3  ifN  =  3 


We  assume  that 

b(t)  e  H  and  t  (t)  e  (L^l'  (Tc)  ^ 

®  ^  0  a  c-  on  Tq 

g  €  L<l(rc) 

where 


(3.6) 

(3.7a) 

(3.7b) 

(3.7c) 


q  =  1  +  max(nin,mT)  ,  q'  =  q/(q-l) 

b(t)  denotes  x  b  (x,t)  ;  t  (t)  denotes  s  -♦  t  (s,t) 

where  s  =  (Sj,  S2, ....  s„)  is  a  point  on  F  (in  panicular  Tp).  Then  the  work  done  by  the 
external  forces  is  defined  by  the  functional, 

F  (t)  €  V' 


<  F  (t).v  )  =  J  b(t)  •  V  dx  +  J  t(t)  •  V  ds 
n  r 

F 

(3.8) 


Also,  to  characterize  the  prescribed  slip  velocity  Uj  on  we  introduce  the  function 


«(t),  6(t)  e  V 

=  U^ft)  a.e.  on  Fq  (3.9) 

For  discussions  of  such  decompositions  of  functions  in  V  into  tangential  and  normal 
components,  see  KIKUCHI  and  ODEN  (12]. 
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Hnally,  to  chancterize  the  work  done  by  nornial  and  firictional  forces  on  Tq,  we  introduce 
the  nonlinear  foims, 

P  :  V  -4  V' 

(  p  (w).  v>  =  f  Cnh(w„)  "  V^ds 
‘C 

j :  V  X  V  ->  m 

Jr 

_  Crh(Wn)'''  Iv^lds 
‘  C 

for  arbitrary  w,  v  in  V. 

We  are  now  ready  to  state  the  following  weak  formulation  of  the  problem. 

Find  the  motion  u(t):  [0,T]  V 
such  that 

<u(t),  v-u(t)>  +a(u(t),  v-u(t)) 

+  C(u(t),  v-u(t))  +  <  P  (u(t)),  V  -  u(t)> 

+  j  (u(t),  V  -  <i>(t) )  -  j  (u(t).  ti(t)  -  <i>(t) ) 

^  <  F(t).  V  -  u(t))  VveV  (3.11) 


(3.10a) 


(3.10b) 


If  djj  (u(t),  u(t)  -  <^t))  denotes  the  partial  subdifferential  of  the  friction  functional  j  with 
respect  to  the  second  argument  (i.e.,  the  velocity),  then  (3. 10)  can  be  written  in  the  equivalent 
form, 

Findu(t)  :  [0,T]  V  such  that 
F(t)  e  (u(t).  (i(t)  -  <I>(t)) 

+  P(u(t))  +  C  u  (t)  +  K  u(t)  +  u’(t)  (3.12) 
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where  K,C  e  o^(V,VO  are  the  operators  defined  by  <  K  w,  v>  =  a  (w,v),  (  C  w,v>  =  c(w,v). 

We  are  now  in  a  position  to  state  a  major  result  on  the  existence  and  uniqueness  of  solutions 
to  (3.11). 

Theorem  3.1  (MARTINS  and  ODEN  [10]).  Let  meas  (r^)>0,  (3.2).  (3.5).  and  (3.6)  hold 
and  suppose  that 

be  l2(0.T;V) 

t  €  l2  (O.T.  ^^(rF))^)  (3.13) 

(so  that  F{t)  6  l2  (0.T;V')) 

Uq€  V,  Uj  e  H 

Then  there  exists  a  unique  solution  u  to  problem  (3.11)  (or.  equivalently.  (3.12))  such  that 

ue  L~(0.r;V) 

u  €  L-  (O.T;H)  n  fO,T;V)  (3.14) 

ii'e  L*(0.T;V') . 

This  result  is  obtained  using  a  regularization  procedure  to  smooth  the  frictionai  boundary  condition 
along  with  a  standard  Galerkin  technique.  It  relies  heavily  on  the  compactness  of  the  nace  operator 
from  V  into  appropriate  boundary  spaces.  The  condition  (3.5),  which  is  used  in  insuring  the 
appropriate  compactness  properties,  is  remarkably  confirmed  by  physical  experiments  on  dry  rough 
surfaces;  see  [9]. 

Let  us  introduce  the  Hilben  space 

w  »{w  «  l2  (0,T;V)  .we  O  (0,T;V’)  ) 
with  narm 

II  W  II  W  »  (  F  ( II  w  (t)  l|2  +  ( II  w  (t)  l|2 )  dt  )*^ 

0  V  H 

An  examination  of  the  proof  of  Theorem  3.1  in  ( 10]  reveals  the  following  continuity  result 


Theorem  3J2.  Let  meas  (Tj^  >  0.  (3.2),  (3.5),  and  (3.6)  hold  and  suppose  that  is  a 

sequence  in  L^(0,T;VO  such  that 

Fk  ->  F  weakly  in  (0,7;^) . 

Then  cortesponding  sequence  if  solutions  u  (Fj^)  of  (3.12)  have  the  property  that 

u  (Fj^)  ->  u  (F)  weakly  in  w 
and  thus 

u  (F|^)  -» u  (F)  in  l2  (0,T;H)  . 

Theorems  3.1  and  3.2  provide  the  basis  for  a  study  of  the  optimal  control  of  distributed 
systems  of  the  type  (3. 1 1).  We  begin  by  introducing  the  notation 

z  e  (0,T;H)  a  prescribed  target  notation  for  t  e  (0,T),  (3. 15a) 

F(t)  e  V 

with 


F(t)  =  (t)  0  .  P  €  H'(O.T)  (3.15b) 

where  0  is  a  prescribed  functional  in  V  defined  by  normalized  external  forces 

<  0,  V  >  *  Jq  b  '  V  dx  Ip  t  •  V  ds  (3. 15c) 

F 

Here  P  s  P(t)  serves  as  a  control  parameter  and  the  functional 


J  :  H*(0.T)  R 


pipinwmyHPfgmw  w  luiit’,  wmjia  wwJWMiuiw'iwwqwipiwKww  ^.-rv’  « 


given  by 

J(P)»(F(llu(-.uP)-z(t)l|2  dt  +  Y  IIPII^  (3.16a) 

0  H  ITCO.T) 

where  u  ( ‘ ,  u  P)  is  the  solution  of  (3.12)  for  the  forcing  function  given  in  (3.16). 

The  optimal  control  problem  that  we  study  here  is  given  as  follows. 

Given  ze  L^(0,T;H)  findPQe  Er(0,T)  (3.16b) 

such  that 

J(Po)  =  infinum  {  J(3) :  p  e  H*  (O.T)  ). 

The  following  existence  result  is  a  consequence  of  Theorem  3.2. 

Theorem  3.3.  Let  (rj))>0,  (3.2),  (3.5),  (3.6)  and  (3.16)  hold.  Then  there  exists  a  solution  to 
problem  (3.16). 

Proof.  Let  (Pil^j  be  a  minimizing  sequence  for  (3.16).  Then  there  exists  a  subsequence  and 
{pj}**  such  that 

Pi  Po  weakly  in  H'(0,T) . 

From  Theorem  3.2  we  see  that 

u  (Pj)  -»  u  (Pq)  weakly  in  l2(0,T;H)  . 

That  Pq  is  a  solution  of  (3.16)  follows  from  the  weak  lower  semicontinuity  of  the  mapping 
P  II  P  IIh’(o,T)  • 

We  have  established  the  existence  of  optimal  controls  for  a  class  of  control  problems 
governed  by  variational  inequalities  modeling  contract  prlblems  with  frictional  effects.  Our  goal 
now  is  to  detenmne  an  approximation  theory  for  these  problems.  Thus,  we  formulate  a  class  of 
approximating  regularized  problems  in  which  the  frictional  condition  is  regularized  similar  to  the 
approach  for  existence  [10].  However,  in  addition  we  include  a  Sobolev  regularization  term.  We 
determine  convergence  behavior  of  the  optimal  controls  for  the  inclusion  of  the  Sobolev  term.  This 
enables  us  to  obtain  reg  irity  properties  for  the  optimal  controls  for  regularized  problems.  These 
properties  allow  us  to  determine  limiting  behavior  for  optimal  control  problems  over  finite 
dimensional  subspaces  of  H'(0,T). 


4.  Regularized  Problems 


In  this  section  we  consider  the  formulation  of  control  problems  governed  by  variational 
problems  that  axe  legulazizations  of  (3.1 1),  c.f.  [10].  In  contrast  to  the  regularization  used  in  the 
papaer  by  Oden  and  Martins  to  establish  the  existence  of  a  solution  to  (3.11),  we  regularize  the 
frictional  contribution  with  a  smoother  function  and  the  introduce  a  Sobolev  regularization  as  well. 
This  is  done  to  permit  us  to  determine  optimality  conditions  for  the  optimal  controls  of  the 
regularized  problems.  Rom  these  conditions  we  may  determine  regularity  properties  of  the  controls 
for  the  regularized  problem.  We  also  obtain  convergence  behavior  of  the  optimal  controls  as  the 
regularizaiton  is  allowed  to  vanish.  We  begin  by  introducing  the  functions 


•Fg:  RN-+R,e>o 
that  satisfy  the  following  properties 


'Fg  €  (R*^,  R)  for  every  e  >  o 

(4.1)(i) 

0  S  'Fg  (v)  S  for  every  £  >  0  and  v  €  R^ 

(4.1)(ii) 

'Fg  (0w  +  (1 . 0)  V)  ^  0  'Fg  (w)  +  (1  -  q)  'Fg  (v) 

(4.1)(iii) 

for  every  e  >  0.  (w,v)  €  x  R^,  and  0  e  [0,1] 

There  exists  positive  constant  D|  and  a  positive  function  e-*  D(e) 

such  that  for  every  e  >  0  and  (w,  v,  u)  €  R^  x  R^  x  R*^, 

(v)  1  ^  Dj  1 V  1 

l'F"g  (w)(v.  w)  1  ^  D,  (e)  1 V  1  (  u  1 

(4.1)(iv) 

There  exists  a  positive  constant  >  0  sucht  that  for  every  e  >  0  and  v  €  R*^ 


I^g  M-lvIl^Dje 


We  define  die  regularized  functional 


jg  :  V  X  V  ->  R 


je  (w.v)  =  Jj-  Cj.  [  (w„  -  g)^f^  'Fg  (v^.)  ds 


We  observe  that  the  derivative  with  respect  to  the  second  argument  of  j,  ( ' ,  • )  is  given  by 


<  Jc  (w,  v),  z  >  =  <  L  (w,  V).  z  > 


so  that 


<  Jg  (w.v),  Z  >  =  Ip  C-r  ( (w„  -  g)J  (v^)  (Zj)ds 


The  variational  formulation  of  the  regularized  problem  of  interest  in  this  study  is  given  as  the 
following. 

Find  a  function  t  —*  Ug(t)  of  [0,  T]  into  V  such  that  for  any  v  e  V 


<  Ug  (t),  v  >  +  e  a  (Ug  (t),  V  )  +  c  (Ug  (t),  v )  +  a  (tig  (t),  v  )  +  <  P  (  Ug  (t),  v  > 
+  <  J.  (u-  (t).  ii,  (t)  -  6  (t) ).  v  >  =  <  f  (i),  v  > 


(4.3)(i) 


with  initial  conditions 


Ug  (0)  =  UQ  6  V 


Ug  (0)  a  Uj  €  V 


(4.3)(ii) 


v"jvvrwjiri.  itj 


•I 


t 

i¥i! 


I;**? 

■‘S 

’k 


’,:i 

‘-!*i 


I 


•::3 


f€  L2(0,  T;  V) 


The  proof  of  the  existence  of  a  unique  solution  to  (4.3)  is  essentially  the  same  as  that  given  in 
[10].  We  use  the  following  convergence  result  the  validity  of  which  follows  firom  a  study  of  the 

proof  of  convergence  of  Ug  to  u  in  [10].  It  should  be  noted  that  for  fe  B  a  bounded  subset  of 
l2  (0,  T;  V’)  that  the  sets 

(Ug  B,e>0} 


{  Ug  (0  f  6  B  ,  e  >  0  } 

are  bounded  in  L"*  (0,  T;  V)  and  L”  (0.  T;  H)  n  (0.  T;  V),  respectively. 

Proposition  4.1.  Let  meas  (fj^)  >  0,  (3.2),  (3.5),  (3.6)  and  (4.1)  hold  and  suppose  that  f, 
-»  f  weakly  in  (0,  T;  V)  and  0  as  i  ^  Then 

Ugj  (fj)  ->  u  (f)  weak  stasji^in  L"  (0,  T;  V) 

Ugj  (fj)  u  (0  weak  star in  L~  (0,  T;  H)  and  weakly  in  (0,  T;  V) . 

Along  similar  lines  but  with  e  >  0  fixed  we  have  the  continuity  result  similar  to  Theorem 
3.2. 

Proposition  4.2.  Let  e>0  be  fixed  and  let  meas  (Fp)  >  0,  (3.2),  (3.5),  (3.6)  and  (4.1) 


hold.  If  ^  f  weakly  in  (0,  T;  VO,  then  (fy)  ->  (f)  weakly  in  W  as  n  ->  oo . 

The  control  problem  for  these  regularized  problems  we  give  as  follows. 

Find  Pj  €  ITCO,  T)  such  that 
J*  (Pp  » infinum  {  J*  (P) :  p  €  ITCO,  T)  ) 

where 

(p)  -F II  Uj  ( • .  t;  p)  -  z  (t)  l|2dt  +  yll  p  l|2 . 

0  H  H’(0.T) 

As  a  consequence  of  Proposition  4.2,  we  have  the  result. 

Corollary  4.1.  Let  the  assumption  of  Proposition  4.2  hold.  There  exists  a  solution  P^  to 
problem  (4.4). 

Remark  4.1.  For  P  fixed  we  note  that 
J'  (P)^J'(P^). 


(4.4) {a) 

(4.4) (b) 


Furthermore,  since  from  Proposition  4.2  we  see  that  Ug  (p)  — »  u  (p)  in  (0,  T;  H),  it  follows 
that 


J*(P)-»J(P) 

as  e  ->  0 .  The  set 

B-  {Pgt  e>0) 


of  optimal  controls  of  (4.4)  is  bounded  in  (0.  T)  and  thus  is  weakly  precompact  in  (0,  T). 


Theorem  4  Ut  meas  (rj5)>0,  (3.2),  (3.5),  (3.6)  and  (4.1)  hold  and  let  Ej 0  as  i 

oo .  Then  any  weak  (0,  T)  limit  point  Pq  of  the  sequence  { is  a  solution  of  (3. 16). 

isl 

Proof.  By  the  above  remaik  there  exists  a  subsequence  again  { Pg^}*  such  that 

ial 

P^->Po  weakly  in  TO.T) 
as  i  — >  <» .  From  Proposition  4. 1  it  follows  that 

l2(0,T;H). 

Thus,  we  see  that 

lim  J®  (Pq)  >  lim  J«'  (p^p  >  J  (P^ 

as  i  oo  and  Pq  is  a  solution  of  (3.16). 

5.  Optimality  Conditions  for  Regularized  Problems 

Having  established  existence  of  optimal  controls  for  regularized  problems  and  their 

convergence  properties  as  e  ->  0,  we  now  seek  to  determine  their  regularity.  We  begin  by  stating  a 
simple  result  that  is  a  consequence  of  the  mean  value  theorem. 


Lemma  5.1.  Let  ^ :  R  — >  R  be  defined  by 
W  » ( (x  -  g)J'" 

for  m  ^  1  and  g  €  R.  Then  for  Xj,  X2  e  R 


■ «.» iji 


1 0  (Xj)  -  ^  (X2)  I  ^  C(in)  ( I  Xj  -  + 1 X2  -  gl™'*)  I  Xj  -  X2 1 . 


Useful  estiinates  for  P  and  j  oiay  now  be  obtained  as  Corollaries  of  Lenuna  5.1. 


Corollary  5.1.  Let  satisfy  (3.6).  Then 


l<P(Ui)-P(U2).v>I^C(m„)[IIUiNliv  +IIu2nIIv  +"5"  „  1 


Proof.  From  Lemma  5.1  it  follows  that 


I  <  P  (u^)  -  P  (U2).  V  >  I  =  I  Ip  Cn  ( [(UjN  -  g)J  "  [(UjN  •  g)J  " )  vn  ds ) 


S  lrC(mn)(luiN*gl  "+lu2N-g)  "  )  lUjN  -  U2N  I  v,^  I  ds 
c 


nin+l  m  -l/m  +1  m  +1  m  -1/m  +1 

S  C(m„)[(JrluiN-gl  ds)  +(lrlu2N-g)  ds)  ] 


.  m-+l  l/m-+l  .  1/m  +1 

(JplUiN'gl  ds)  (Iplv^l  ds)  ] 


i| 

I’l 

y 

y 


I 

I 

il 

p’' 

U 

E*' 

m\ 


\ 


Hence,  we  see  that 


I  <  P  (Uj)  -  P  (U2),  V  >  I  ^  C  (m„)  [  II  u,N 


mn-l  m  -1 

I  +IIu2mII  ] 

L"’n‘^*(rj  L"’n+'(rj 


■llfl*  s 

+  II  g  II  ]  IIUlKI  -  U^K,  II  llv^l 

L  "*n+*  (Fg)  ^  L  '"n**  (r^)  L  (F,.) 


Recalling  that,  under  the  assumption  (3.6),  imbeds  continuously  into  L'l(r  )  and 

c 

V«H*(Q)  nu^)S  continuously  into  we  have  the  result 

In  a  similar  manner,  we  have  for  j. 

Corollary  5J.  Let  mj.  satisfy  (3.6).  Then 

.  mT*!  tnyl  fn-T^l 

Ij  (Ui.v)  -  j  (Uj,  V)  i  ^  C  (nvr)  [  Iluj  lly  +  IIU2  Hy  +  II  g  II  ]  Hu,  -  uJIv  llvllv 

L  (Tg) 

Let  h  €  L^  (0,T)  and  define 

wg  =  (u  (3  +  5h)  -  \)  (3) )  /  5  for  5  >  0 . 

We  note  that  w  =  wg  satisfies 

(w  (t).  V)  +  e  a  (w  (t),  V)  +  c  (w  (t).  V)  +  a  (w  (t).  v)  + 

+  1/6  <  P  (u  (t;  3  +  5h)  -  P  (u  (t;  3)).  v  >  +  (5j) 

+  1/5  (<  J  (u  (t;  3  +  5h).  (u  (t;  3  +  6h)  -  O  (t)).v  >  -  <  J  (u  (t;  3).4^u  (t;  3)  -  6  (i)),v  >  ) 

=  h  (t)  <  0,  V  > 


We  now  use  Corollary  5. 1  to  obtain  for  5  >  0  the  estimate 
1/5 1  <  P  (u  (t;  3  +  5h)  -  P  (u  (t;  3))  -  v  >  I  ^ 

ni„-l  m„-l  m_-l 


n  II  * 

^  C  (in„)  ( llu  (t;  3  +  5h)lly  +  Hu  (t;  3)llv  +  H  g  II  )  llw(t)llv  Hvllw 

L  '"n+‘  ^  ^ 


and 


1/5 1  <  J  (u  (t;  3  +  6h),  iid;  3  +  5h)  -  O  (i)),  v  >  -  <  J  (u  (t;  3),  u  (t;  3)  -  O),  v  >  I 


1/5  Jj.  CfC  [(u„)  (t;  3  +  Si)  -  g) J  <f’(u  (t;  p  +  5h)  -  O  (t)  (v)  - 


[(Un)  3)  -  g)J  (t:  p)  -  6  (t)  (V))  ds 


Now  set 


A  *  g  *  ^  ([(Un)(‘i  P  +  5h)  -  g;^I  -  [^u  (t;  B)  -  g)  J  '^’(u  (t;  P  +  6h)  -  <l>(t))(v)  ds 


B  =  -^1  Ir  Cj.  [(Un)(t;  p)  -  g)^]  (t;  p  +  5h)"^C>(t))  -  OXii  (t;  P)  -  6(t))(v)  ds 


From  (4, 1)  (w)  we  see  that  for  A  we  have 


j>  m-T- 1  m-n*  1 

A  ^  Dj  I J j<  Cj-  mj(  lu„  (t;  p  +  6h)  -  gl  +  |u„  (i;  P)  -  gl  )  Iw  (i)l  (v)  ds  I 


and  for  B 


B  ^  Dj(e)  I  Jp  Cj  lu^  (t;  P)  -  gl  Iwl  (v)  ds 


We  now  make  the  restriction  of  assumption  (3.6) 

l^Snvj.  <  +  «»ifN  =  2 
i  2  if  N=3 

Under  the  assumption  (5.2)  we  have  chat 


A  *  rn*T'‘^2  fn-p/ni*T^2 

B  ^  D|(e)  Cp  ( Jp  lu„(t;  P)  -  gl  ds)  ll\^l  llvll 

c  L  '"n*2(r^)  L 


^  C  (e.  fl»r)  ( liu  (B)llv  +  llgll  )  llwll  V  livil  V 

^  L  “‘n+kr^) 


Thus,  we  see  that  from  (S.l)  -  (S.3)  and  (S.S), 


(w(t),  V)  +  e  a  (w(t),  v)  +  c  (w(t),  v)  +  a  (w(t).  v) 

^  I  h(t)  I  <  e.  V  >  +  C  (5)  ( llw(t)llv  +  llw(t)llv  )  llv|ly 


where  C(5)  depends  on  the  parameter  of  the  problem  as  well  as  u(3),  g,  and  e  and  is  bounded  for 


S>0  in  a  bounded  set  Now  setting  v  s  w(t)  we  have 


d/dt  ( llw(t)l|2  +  e  a  (w(t),  w(t))  +  a  (w(t),  w(t))  +  2C  (w(t),  w(t)) 
H 


^  2 1  h(t)  1  <  e,  w{t)  >  +  2^  (5)  ( llw(t)llv  +  llw(t)llv  +  llw(t)|lv  ) 


^  C,  I  h(t)l  +  Cj  (5)  ( llw(t)llv  +  llw(t)llv ) 


We  may  now  obtain  from  Gromwall's  inequality  that 


Ilw5(t)l|2^+  Ilw5(t)l|2^+  J'  llw5(T)l|2^t  ^  C(5)  lh(T)|2  dt 


where  C(5)  is  a  function  of  5  sending  bounded  sets  in  iR'*'  into  boundet  sets  in  ^[iTa  similar 
manner,  using  v«w(t),  we  may  show  boundedness  of  w(t)  in  l2(0,T;V).  Hence,  we  see  that 


there  is  a  sequence  5.-)0  such  that 


Wj.->  w  weak  star  in  L"(0,T;V) 


star  in  L’"(0,T;V) 


Wj.-*  w  wtakj^r  in  L^(0,T;V) 


u 

With  this  preparatioii,  we  have  the  following  result:  We  denote  by  dug(b)(h)  the  variai^of 
Ug(P)  with  increment  h. 

Proposidon  5.1.  Let  meas  (Fjj)  >0,  (3.2),  (3.5),  (3.6),  (4.1)  and  (5.4)  hold.  Then  w  =  9ug(p) 
(h)  exists  and  satisfies  the  variational  equadon 

(w(t),  V)  +  ea(w(t),  V)  +  C(w(t),  V)  +  a(w(t).  v)  +  /  ((u„  -  g)^  w^d),  v„  ds  + 

+  J  {  Cj  [(u„  -  g)J*^ (Uj  -  6)  (w,.,  Vj)  -»•  (5.10)(a) 

+  m|.  Cj.  [(u„  -  g)J  ^  'F*  (uy  -  <P)  (v^.)  }  ds  =  h(t)  <  ©.  v  > 

for  any  V  €  V  and  with  initial  condition 

w(0)«w(0)«0.  (5.10)(b) 

Based  on  estimates  similar  to  (5.7),  it  is  straightforward  to  obtain  the  following  result 

Lemma  5.2.  Let  (4.1)  hold.  Then  there  exists  a  unique  solution  w  to  problem  (5.9). 

Proof  (of  Proposidon  5.1).  From  (5.8)  it  is  clear  that  the  first  four  terms  of  (5.1)  converge  to  the 
first  four  terms  of  (5.9)(a).  For  the  boundary  terms  we  use  the  compactness  cf  [10]  to  see  that 

Wj .  ^  w  in  L^(0,  T;  L*!  (Tg)) 

and 

wj .  -4  w  in  L^O,  T;  L*!  (f^,)) 

Furtheimoie,  from  (5.7)  it  follows  that 
u(p  +  5h)  -+  u(P)  in  H  (0.  T;  V) 


and  therefore 


uO  +  5h)  ->  u(P)  in  L^CO,  T;  L**  (Tg) 
where  l^ii^forN«3  and  forNs2. 

Having  observed  the  above  convergence  properties,  we  note  that  the  boundary  integral  terms 
convey  by  the  dominated  convergence  theorem. 

Remark  5.1.  We  point  out  that  the  Sobolev  regularization  is  used  to  obtain  inequality  (5.S).  It  is 
necessary  since  the  coefficient  of  the  ItwCOIIy^  term  in  inequality  (5.7)  cannot  be  made  small 
Using  the  results  from  Proposition  S.l.  we  may  now  calculate  the  variation  of  J. 


Corollary  53.  Let  meas  (Fjj)  >  0,  (3.2),  (3.5),  (3.6),  (4.1)  and  (5.4)  hold  and  let  J*  be  given 
by  (4.4)(b).  Then  9J*(p)(h)  is  given  by 

aj«(P)(h)  =2  jTjUg  ( t;  3)  -  z(t),  w(t,h))^  dt  +  2  y(3.  h)H.(o  t)  (5- 1  D 

where  w  is  the  solution  of  (5.10). 

Having  determined  the  differentiability  of  the  functional  we  have  the  result 


Proposition  5.1.  Let  meas  (f^)  >  0,  (3.2),  (3.5),  (3.6),  (4.1)  and  (5.4)  hold.  If  3^  is  a 
solution  of  (4.4),  then 

/T  (Uj  (t;  3)  -  z(t).  w(t,h))^  dt  +  Y(3g,  h)H(o.-n  =  0  (5.12) 


for  all  h  6  l2(0,T). 

We  may  now  use  Proposition  5. 1  to  obtain  regularity  of  3^.  E  >  0 . 

Theorem  5.1.  Let  meas  (Fp)  >  0,  (3.2),  (3.5),  (3.6),  (4.1)  and  (5.4)  hold.  If  3^  »s  a  solution  of 
(4.4),  then  h2(0,T). 


Proof.  Consider  the  linear  mapping 


h  -*  a(b)  ■  (t;  -  z(t), 


of  frCO,  T)  into  R.  We  note  from  Proposition  4.2  and  from  (5.8)  that 

I  a(h)  I  ^  II  Ug  ( • .  pg)  -  z  ( • )  II  l2<0.T;H)  C  M  h  II  j  2(o,t) 

and  thus  in  a  continuous  linear  function  on  L^(0,T).  It  follows  by  the  Riesz(?)  Representation 
theorem  [16]  that  there  exists  a  6  L^(0,T)  such  that 
a(h) » (a.h)  l2(o.t)  • 


From  (5.12)  we  see  that  for  every  h  €  H'(0,T) 

(Pji  h)H’(0,T)  “  *  ■  (5.13) 

We  conclude  from  equation  (5.13)  that  in  fact 

-  P«.  *  Pe  =  •  ■'If® 

with 

Pe(0)»Pe(T)  =  0. 

Hence.  3^  €  h2(0.T). 


6.  Approximation 

We  begin  by  giving  the  Galerkin  formulation  for  problem  (4.3).  Thus,  let  C  V  with 

basis  (  B:  and  let 

i-l 

u^  (x,t)  ■  ®i 


be  such  that  for  e  2  0  fixed  and  v  € 


<  V  >  +  e  a(u^(t),  V)  +  c(Ug  (t),  v)  +  a(Uj  (t),  v)  +  <  P  (u^  (t),  v  )  + 

+  ( (Ug  (t).  Uj  (t)  -  f  (t)).  V  )  =  3(0  <  e.  V  >  (6.  l)(a) 

with  initial  coodititms 


N 

Uj  (0)  *  Uq 

'N N 

Uj  (0)  =  Uj 


(6.1)(b) 


where  Uq  and  U|  converge  to  Uq  and  Uj  in  V,  respectively,  as  N  ->  <» .  The  following  is  not 
difficult  to  show. 


Proposition  6.1.  Let  meas  (fj^)  >  0.  (3.2),  (3.5),  (3.6),  (4.1)  and  (5.4)  hold  and  let  e  >  0  be 

fixed.  Let  13})**  a  sequence  in  H’(0,T)  such  that  3i  P  weakly  in  ^(0,7)  and  let  Nj  «<►. 
ial 

Then 

N, 

Uj  (3i)  Uj  (3)  in  l2(0,T;H) 


as  i  — ♦  o®  . 

To  approximate  the  controls,  let  L^  be  a  subspace  in  H'(0,T).  Suppose  there  is  a 
continuous  linear  mapping  H'(0,T)  -» L^  such  that 

II  u  -  llif(o,T)  ^  C II  u  IIh'(o,t) 
and  for  u  €  H^O.T) 


II  u  -  ••h'(0,T)  ^  “  *’^(0.7) 
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where  and  5(M)^  with  5(N)— »0  as  M—*o» .  Both  C  and  5(M)  are  independent  of  u. 
The  proUems  we  consider  are  given  by  the  following; 

Rnd  such  that 


Jm  0“'")  -  Win™  1 (P) :  p  e  lN) 

where 

Jm  (3)  =  (tip)  -  z(t)  Y II 3  if^Q  ^ } 


(6.3)(a) 


(6.3)(b) 


Remark  6.1.  The  existence  of  solutions  of  (6.3).  Moreover,  these  solutions  are  bounded 

t 


in  H'(0,T). 

Theorem  6.1.  Let  be  a  weak  ITCO.T)  limit  point  of  Pg  '  ‘  for  — >  «>o  and  Nj 

Then  Pg  is  a  solution  of  (4.4). 

Proof.  Since  Pg  -»  pg  weakly  in  H'(0,T)  as  i  ,  it  follows  from  Proposition  6. 1  that 


M:.N:  M:,N: 

“c  (Pe  )->  Ue(Pe) 
in  l2(0,T)  as  i  — >  oo .  Thus,  we  see  that 


N:  M:.N: 

Jg  (Pg  )^Jg(Pg) 

Now  from  Theorem  5.1,  Pg  €  H'(0,T)  and,  thus  from  (6,2),  we  see  that 


Nj  Mi  N:  M:  N: 

J,  ‘a  p,)  2  J,‘  (P,'  “) 


and 


N:  M:  N, 

Jg(pg)  -lim  Jg  (pg  )  >  JgCPg)  . 
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•T*ir 


We  conchide  that  is  t  solution  of  (4.4). 

Remark  <S.l.  It  is  now  clear  that  an  iterated  limit  theorem  follows  from  Theorem  4.2  and  Theorem 

6.1. 

Theorem  62,.  Let  ->  0.  Then  there  exist  sequences  Mj  Nj  -♦(»•.  ««•)  as  i  -♦  o*  such  that 


M^Nj 

w-lim  (w-lim  ,  )=  Po 

ej-*0  T 


where  P^is  a  solution  of  (3.16)  and  where  the  w-limits  are  weak  limits  in  pr(0,T).  Also, 

MiNj 

lim  Urn  Jp.  (p,  )  =  J  (Pq)  . 


Acknowledgement.  The  suppon  of  this  work  by  the  U.  S.  Air  Force  Office  of  Scientific  Research 
under  Contract  F49620-86-<i-0051  under  grant  No.  AFOSR-84-027 1  is  gratefully  appreciated. 


REFERENCES 


1.  Duvaut,  G.  and  Lions,  J.L.,  Inequalities  in  Mechanics  and  Physics,  Springer- Verlag, 
Beilin,  Heidelberg,  New  York,  1976. 

2.  Ldtstedt,  P.,  "Coulomb  Friction  in  Two-Dimensional  Rigid  Body  Systems",  ZAMM,  61,  pp. 
605-615,  1984. 

3.  Jean,  M.,  and  Pratt,  E.,  "A  System  of  Rigid  Bodies  with  Dry  Friction",  International  Journal 
of  Engineering  Sciences,  Vol.  23,  No.  5,  pp.  497-513,  1985. 

4.  Schatzmann,  M.,  "A  Class  of  Nonlinear  Differential  Equations  of  Second  Order  in  Time, 
Nonlinear  Analysis,  Theory,  Methods  and  Applications,  2,  pp.  355-373, 1978. 


26 


I 


5.  Canoo,  M  and  Pascali,  E.,  ”11  Pioblema  del  Rimbalzo  Unidimensionale  e  sue 
A^^nossiinazioni  c<»  Penalizzaziom  non  Convesse",  Rend.  Mat.,  (6)  13,  4,  pp.  541-553, 

6.  Ldtstedt,  P.,  "Mechanical  Systems  of  Rigid  Bodies  Subject  to  Unilateral  Constraints",  SIAM 
Journal  of  implied  Mathematics,  Vol.  42,  No.  2,  pp.  281-296,  1982. 

7.  Amcrio,  L.  and  Prouse,  G.,  "Study  of  the  Motion  of  a  String  Vibrating  Against  an  Obstacle", 
Rend.  Mat.,  (2)  8.  Ser.  VI,  pp.  563-585,  1975. 

8.  Ameiio,  L.,  "Su  un  Pioblema  di  Vincoli  Unilaterali  per  lEquazione  non  Omogenea  della  Corda 
Vibiante,  lAC  (Inst  per  le  Applicazioni  del  Calcolo  "Mauro  Picone")  Publicazioni,  Scr.  D 
109,  pp.  3-11,  1976. 

9.  Oden,  J.  T.  and  Martins,  J.A.C.,  "Models  and  Computational  Methods  for  Dynamic  Friction 
Phenomena",  Computer  Methods  in  Applied  Mechanics  and  Engineering,  52,  pp.  527-634, 
1985. 

10.  Martins,  J.A.C.  and  Oden,  J.T.,  "Existence  and  Uniqueness  Results  for  Dynamic  Contact 
Problems  with  Nonlinear  Normal  and  Friction  Interface  Laws",  Journal  of  Nonlinear  Analysis 
(to  appear.) 

11.  Rabier,  P.  J.,  Manins,  J.A.C.,  Oden,  J.T.,  and  Campos,  L.,  "Existence  and  Local 
Uniqueness  of  Solutions  to  Contact  Problems  in  Elasticity  with  Nonlinear  Friction  Laws", 
International  Journal  of  Engineering  Sciences,  Vol.  24,  No.  11,  pp.  1755-1768,  1986. 

12.  Kikuchi,  N.  and  Oden,  J.T.,  Contact  Problems  in  Elasticity,  SLAM  Publications, 
Philadelphia,  1987. 


13.  Back,  N.,  Burdekin,  M.  and  Cowley,  A.,  "Review  of  the  Research  on  Fixed  and  Sliding 
Joints",  Proc.  13th  International  Machine  Tool  Design  and  Research  Conference,  ed.  by  S.A. 
Tobias  and  F.  Koenigsberger,  MacMillan,  London,  pp.  87-97,  1973. 

14.  Back,  N.,  Burdekin,  M.  and  Cowley,  A.,  "Analysis  of  Machine  Tool  Joints  by  the  Finite 
Element  Method,"  Proc.  14th  International  Machine  Tool  Design  and  Research  Conference, 
ed.  by  S.  A.  Tobias  and  F.  Koenigsberger,  MacMillan,  London,  pp.  529-537,  1974. 

15.  Lions,  J.L.,  Quelques  M4thodes  de  Resolution  des  Problimes  aux  Limites  Non 
Lindaircs,  Dunod,  Paris,  1969. 

16.  Luenberger,  D.G.,  Optimization  by  Vector  Space  Methods,  Wiley,  New  York,  1969. 


27 


T^-TJ 


Appendix  B 


Regularity  properties  and  strict  complementarity 
of  the  output-least-squares  approach  to  parameter 
estimation  in  parabolic  equations. 
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1.  Introduction. 


In  this  paper  we  study  properties  of  the  output-least-squares 
formulation  for  the  estimation  of  the  diffusion  coefficient  a  in 

(1.1)  Ut  =  (au,),  +  f, 

together  with  appropriate  initial-  and  boundary  conditions,  from  data  z. 
More  precisely,  assume  that  (1.1)  is  a  model  of  a  system  S  where  the 
coefficient  a  is  unknown.  To  determine  a  we  assume  the  availability  of 
data  z  of  S.  The  task  arises  of  determining  a*  such  that  u(a*)  best 
fits  the  data  z.  The  moat  comoion  approach  to  this  parameter  estimation 
problem  is  the  output  least  squares  formulation,  given  by 

(P)  ainiaize  lu(a}  -  zi^ 

over  some  set  of  admissible  parameters  A,  where  u(a)  satisfies  (1.1)  and 
the  distance  u(a)  -  z  is  taken  in  an  appropriate  norm.  This  problem  has 
received  a  considerable  amount  of  attention  (see  [3,4,5,11,12,18]  et  a^. 
These  investigations  are  primarily  concerned  with  the  approximation  of  the 
infinite  dimensional  problem  (P),  with  numerical  aspects  of  the  finite 
dimensional  approximating  problems  as  well  as  with  establishing  techniques 
to  cope  with  the  illposedness  of  the  estimation  problem,  by  which  we  mean 
the  lack  of  continuous  dependence  of  the  solutions  of  (P)  on  perturbations 
in  z.  Predominantly  among  the  latter  is  the  regularization  technique,  in 
which  (P)  is  replaced  by 

(P)®  ainiaize  lu(a)  -  zi*  Blal*,  over  A, 

where  S  is  a  small  positive  parameter  and  I  ■  I  denotes  an  appropriate 
norm.  It  appears  that  optimization  theoretic  aspects  of  (P)  and  (P)®  have 
not  yet  been  studied  in  detail  and  we  shall  make  a  step  in  this  direction. 

Two  constraints  will  be  involved  in  defining  A:  a  pointwise  lower 
bound  on  the  elements  in  A  guarantees  that  the  differential  operator 
appearing  in  (1.1)  is  elliptic  and  a  norm  bound  is  used  to  establish 
existence  of  solutions  of  (P).  Due  to  the  lack  of  radial  unboundedness  of 
the  oaapping  a  ■*  u(a),  (P)  may  have  no  solution  unless  the  elements  in  A 
are  bounded  in  an  appropriate  sense  or  unless  a  regularization  term,  as  in 
(P)®,  is  used. 


In  Motion  2  we  collect  reeulta  on  reffularity  propertiee  of  the  solution 
of  (l.l)*  Special  attention  is  fiven  to  the  fact  that  theM  results  are 
obtained  under  weak  smoothness  requirements  on  the  coefficient  a.  In 
sections  3  and  4  we  show  that  due  to  the  norm  constraint  or  the 
regularization  term,  the  solutions  a  of  (P)  experience  certain  local  as  well 
as  global  regularity  propertiest  by  which  we  mean  that  a  is  more  regular 
than  the  functions  in  the  set  of  admissible  parameters  A*  While  section  J  is 
devoted  to  the  csm  where  a  is  a  function  of  the  spatial  variable  only,  we 
treat  space  and  time  dependent  variables  in  section  4.  Besides  being  of 
interest  in  their  cwn  right,  regularity  properties  of  this  kind  proved  to  the 
uMful  for  elliptic  estimation  problems:  in  [12]  rate  of  convergence  results 
for  Galerkin  approximations  of  the  output-least-squares  formulation  of 
estimation  problems  was  proved  and  this  rate  depends  on  the  regularity  of 
the  solutions  of  the  infinite  dimensional  problem,  and  in  [14]  a  priori 
knowledge  of  the  regularity  of  the  solution  improved  results  for  the 
penality  oieihod  formulation  of  the  norm  constraint. 

In  section  5  we  give  conditions  under  which  the  norm  constraint  is  in 
fact  active  and  then  we  produce  a  claM  of  examples  for  which  the 
formulation  of  (P)  without  norm  constraint  on  A  has  no  solution.  The 
proof  of  the  activeness  of  the  norm  constraint  is  achieved  by  showing  that 
the  asaociated  Lagrange  multiplier  is  nontrivial.  Prom  comparision  with 
estiaiation  problems  for  elliptic  equations  we  expect  that  nontriviality  of  the 
Lagrange  multiplier  will  have  important  conMquences  for  parabolic 
equations  as  welL  Among  them  are  the  stability  of  the  solutions  of  (P)  with 
respect  to  perturbations  in  z.  This  will  be  studied  elMwhere. 

Throughout  we  um  standard  notation.  The  norms  of  inner  products  of 
elements  in  function  spaces  are  indicated  with  a  subcript,  as  for  example 
I'lgj.  However,  if  the  function  space  is  L*  then  we  nay  drop  its  notation 
and  we  use  <  * ,  *  >  for  the  inner  product  in  L*  and  I  ■  I  for  the  norm 
in  L*  and  in  R". 


2. 


In  this  section  we  summarize  some  results  concerning  the  output  least 
squares  and  the  regularized  output  least  squares  formulation  of 


(2.1) 


=  (au,),  f  in  Q 
u(x,0)  =  ♦(x)  in  Q 
u(O.t)  =  u(l,t)  =0  for  t  c  (O.T) 


where  Q  =  (0|1)  and  Q  =  Q  a  (0,T),  T  ^  0.  Depending  on  the  cases  when 
a  is  a  function  of  x  only,  or  a  function  of  x  and  t,  we  define  the 
admissible  parameter  sets 


^  =  (a  =  a(x):  a  e  H*(0,1),  a(x)  fc  i/  0,  laljji  *  n) ,  v  ^  ii  *  • 
or 

^  =  (a  =  a(x,t);a  €  H*(0,1)  ,a(x,t)  a  i/  ^  Oilal^,  *  n) ,  vT  ^  *  • 

The  requirements  0  ^  v  ^  n  and  Tv  ^  u  guarantee  that  A  is  not 
empty.  The  set  A  is  given  the  topology  of  H‘  or  H*  respectively.  We 
shall  employ  the  following  fit-to-data  criteria: 


and 


Ji(a) 


T  I 

=  J  J  lu(x,t;a)  -  Zi(x,t)  I  *dxdt  + 

O  0 


T 

Ja(a)  =  J  lu(x,T;a)  -  2a(x)l*dx  +  B'al*, 


where  B  a  0,  Zi  c  L* (Q) ,  Za  c  L* (Q)  and  I  ■  l/|  stands  for  *  ’ ' gi  if  we 
consider  only  x-dependent  coefficients  and  for  ' '  *  ga  if  the  coefficients 
depend  on  x  and  t.  Finally  u  =  u(x,t;a)  denotes  the  solution  of  (2.1) 
depending  as  a  fuention  of  x  and  t  on  the  coefficient  a.  The  least 
squares  problems  that  we  shall  investigate  are  given  by 

(Pf)  ain  J4(a)  over  A. 

By  putting  it  s  •  or  0  =  0  this  formulation  includes  the  case  without 
norm  bound  on  the  parameter  set  or  without  regulariaation  in  the  cost 
functional.  If  the  dependence  of  J|  or  (P|)  on  B  or  it  is  relevant 


uunuuKiitiivmfliuvuvvii 


JTWWJTn 


i  ITS  TTil  irm 


we  denote  this  bj  additional  indices,  as  for  instance  (P,  )A.  First  some 
existence  and  regularity  properties  for  the  equation  (2.1)  are  summarized. 
Unless  specified  otherwise  it  is  assumed  that  #  e  L*(Q),  f  e  L‘(0,T;ir*) 
and  a  «  A.  Henceforth  we  shall  write  L’(X)  for  L’(0,T;X)  when  X  is  a 
Banach  space.  Let 

W(0,T)  =  (u:  u  €  L*(HM,  u*  e  L*(ir‘)) 

o 

endowed  with  *  *“**l*(H“‘)^*^*  **  *  nor*.  It  is  wel^mown  that 

W(0,T)  is  a  Hilbert  space  which  embeds  continuously  into  C(L*)  [19,p.  19] 
and  compactly  into  L*(C).  Here  we  used  the  following 


i.^iiime  g.i  [22,  p.  271].  Let  Xo,  X,  X|  be  Banach  spaces  such  that  Xo  ^  X 
c  X|,  where  the  injections  are  continuous,  X|  are  reflexive  for  i  =  0,1 
and  X«  embeds  compactly  in  X.  Set 

Y  =  (v  €  L*(Xo)  :  v’  €  LMX,)), 

with  Ivl*  s  '''’l*(X  )  ^  ''^**L*(X  )'  injection  of  Y  into  L*(X) 

is  compact. 


Proposition  2.1.  For  a  =  a(x,t)  c  A  there  exists  a  unique  solution  u  of 
(2.1)  in  W(0,T). 


Proof.  Existence  end  uniqueness  of  u  c  L*(H‘)  is  shown  in  [16,  p.  121]. 
Usinf  this  in  (2.1)  we  find  u  <  W(0,T).  lioreover  we  have  the  estunate 

(16,  p.  118]. 

(2.2) 

a 

where  k  depends  on  v  and  p  but  is  otherwise  independent  of  a. 

ProposiUon  2.2.  If  #  e  L*(Q),  f  c  L^L*)  and  a  is  time  independent  so 
that  a  :  a(x)  c  H'(0),  s(x)  *  v  ^  0,  then  the  solution  of  (2.1)  satisfies 
u  c  L*(C*). 


Proof.  Uaing  a  aemigroup  theoretic  argument  let  A  in  L*  be  defined  by 

dom(A)  =  H*  and  Av  =  (av.),.  It  is  twllknown  that  A  is  the 

o 

infinitesimal  generator  of  an  analytic  semigroup  S(t)  and  that  the  eolution 
u  of  (2.1)  can  be  expreaaed  as 

t 

(2.3)  u(t)  =  S(t)#  +  J  S(t-s)f(s)ds. 

o 

In  the  following  estimates  we  use  the  fact  that  the  fractional  power  space 
doa(A‘*)  endowed  with  lA^vl  as  a  norm  embeds  continuously  into  C*  if 
a  e  (3/«tl)i  C9iP.  39].  Recall  that  IA‘*S(t)l  *  t~<’‘  for  f  ^  0.  Since 

T 

IS(t)#l^i^^,^  ^  ^  1  'A<*3(t)4ldt  ^  •,  where  k  is  an  embedding  constant, 

we  find  that  S(t)#  e  L‘(CM.  Moreover  using  the  convolution  theorem  we 
have 

T  t  T  t 

J  J  IS(t-a)f(s)t^idsdt  *  k  f  f  l A«S(t-s)f(s) Idsdt 
«  0  0  0 

T  t 

s  J  J  s~®dsdt 

'  '  o  o 

4  •  ^  -  If  I  Ta-a 

Using  these  estimates  in  (2.3)  implies  the  claim. 


Proposition  2.3.  If  f  e  L*(L*),  ♦  «  HI  and  a  =  a(x,t)  e  A  «  H*,  then 


u  c  L*(H*  n  H*),  Ut  €  L*(L*),  and  u  e  C(H»). 

o  o 


Proof.  This  result  is  stated  without  proof  in  [17,  p.  180].  Since  the  other 
wellknown  references  require  stronger  (W*>*  -)  regularity  assumptions  on 
a  and  since  regularity  of  a  is  the  focus  of  our  considerations,  we  outline 
the  proof. 

Employing  the  mean  value  theorem  for  integrals  one  can  prove  that 


I 

(2.4)  j  ^  *'*'h*(0) 

0 


and 


t 

(2.5)  (I  la,(x,  )l*.(g  .jjdx)’/*  -  *'*'h»(Q)' 

o 

for  a  conatant  K  independent  of  a  e  H'. 

We  recalled  before  that  u  e  L*(H*).  Using  (2.4)  we  can  now  refer  to  [17, 

0 

p.  178]  to  conclude  that  u  e  L*(H‘)  imd  e  L*(L*).  Next  observe  that 
observe  that  (au„)y  -  a„u,  =  -  f  -  a„u,  and  the  right  hand  side  of 

this  expreaaaion  is  in  L^(L‘).  Consequently  au„x  =  (suy),  -  SkU,  is  in 
L*(L‘)  and  u  e  L*(W*»‘).  Recall  that  W»»‘  embeds  continuously  into  C. 
This  implies  that  u  e  L*(C‘)-  Next  an  estimate  for  u*,  in  L*(L*)  can  be 
derived; 


*  '“^'l*(L*) 


T  » 

*  J  la,u,l  »dxdt) 

0  o 

T  I 

+  (J  ''**^  '^^'l*(0,l)  J  ^ 'l*(0,T)*"**^*^* 

0  9 

+  K'u'La(ci)'®'ga. 


where  we  used  (2.5)  in  the  last  step.  Thus  u  e  L'(H*).  Employing  this  fact 
in  (2.1)  we  obtain  Ut  c  L’(L’).  The  final  claim  follows  from  [19»  p.  19]. 


To  develop  the  final  regularity  result  that  will  be  required  we  need  to 
express  (2.1)  in  an  evolution  operator  theoretic  setting.  For  t  e  [0,T]  let 
the  operators  A(t)  be  given  by 


dam  A(t)  =  H»  «  H‘ 


with 


A(t)#  =  (a(t,  )♦,),  for  ♦  e  H»  n  H‘, 
and  for  A  c  C  let 


R(X;A(t))  =  (A  -  A(t))-', 


whenever  R(X;AU))  c  £(L*(0)).  Here  £(L’(Q))  denotee  the  set  of  bounded 
linear  operatori  on  L*(Q).  The  following  lemma  will  be  required. 


Lemma  2.2. 

(a)  If  a  €  then  a  e  C‘ (Q)),  which  is  the  space  of  Holder 

continuous  functions  with  exponent  1/2  and  values  in  .  Let  a  =  a(x,l)  e 

A. 

(b)  There  exist  constants  c  and  c  such  that 

-  '♦'h^CQ)  *  *'  ^  '^'h*(Q)’ 

for  all  ♦  e  n  H‘  and  t  e  tO,T], 

(c)  ^  7^  for  all  A  e  C  with  Re  X  k  0  and  t  e  C0,T], 

(d)  there  exists  a  constant  L  such  that 

l(A(t)  -  A(s))A(t)-‘I^^j^,,qjj  *  LIt-sl*/* 
for  all  t,  s,  re  [O.Tl. 

From  the  proof  it  can  be  seen  that  c,  c  and  L  in  (b)  and  (d)  can  be 
taken  independentl7  of  a  e  A»  if  p  in  the  definition  of  A  is  finite. 


Proof  of  Lemma  2.2. 

(a)  Let  t,  s  e  (0,T].  Then 

t  >  t 

[  J  laH(t,x)  -  a*(s,x)l*dx]‘^*  «  (  JlJ  a,t(fftx)dal *dx] 

0  0  • 

1  1 

s  lt-s»‘/*  [  J  J  »a,t(»tX)l*dsdx) 

0  0 

and  thus  (a)  follows. 

(b)  The  second  inequalit7  in  (b)  is  obvious.  The  first  estimate  is  fairl7 
standard  as  wall  and  we  onl7  sketch  the  proof  to  convince  the  reader  that 
S  can  be  taken  uniforml7  with  respect  to  t  c  [0|T].  Suppressing  the 
dependence  of  s  on  x  we  find 


iUimumWMUllUi-'m  j  M  wwum  .■mry  ■■>  .■^Ii;wraw;fwy\fwwyv 
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Ill  I 

J  8*(t)#*,dx  <  I  (A(t)#)*dx  -  2  J  a(t)a«(t)#,#xK<lJ‘  “  J  a*#Kdx 

0  0  0  0 

1  1  t  1 

*  f  (A(t)#)*dx  +  5  J  a*(t)#xxdx  +  2  J  ax(t)#Jdx  -  J  ax(t)#xdx. 

0  O  0  0 

Hence  it  followa  that 

I  1  1 

~  j  a*(t)#xxdx  <  I  (A{t)#)*dx  +  J  ax(t)#xdx 

0  0  0 

and  ainca  a  €  A 

1  1  1 

1/  J  ♦*xdx  *  2  f  (A(t)#)»dx  +  2l#xl*.  J  ax(t)dx. 

0  0  o 

From  (a)  and  the  Ehrlingache  Lemma  [23,  p.  118]  there  exiata  for  every 
e  ^  0  a  conatant  c(e)  auch  that 

a  2lA(t)^l»  +  2(€I#x,I»  +  c(e)l#,l*)lalg,^^j. 

For  an  appropriate  choice  of  c  thia  impliea 

(2.6)  •♦xJ  *  clA(t)#l 

for  all  4  e  H*  n  H*  and  t  e  [O.T].  The  eatiaate  '^Iqi  *  c(A(t)#l  for 

aoma  conatant  c  uniforaly  in  f  c  H*  n  and  t  e  [0,T]  ia  eaaily 

0 

verified  and  together  with  (2.6)  thia  iapliea  (b). 

(c)  For  #  c  H*  n  we  have 

O 

(♦,A(t)#)  a  -vlgx»», 

and  thus  aemigroup  theory  impliea  that  for  every  t  e  [0,T],  A(t)  ganeraiea 
a  semigroup  exp  A(t)  on  L*(0)  aatisfying  lexp  A(t)l  a  a'*^^  [20;  pp. 

8,14].  This  implies  (c)  (compare  [20,  p.  11]). 

(d)  The  operator  A(r),  r  e  [0,T]  are  homeomorphiaas  bett«een  H*  n  H*  ond 
L*  and  in  view  of  (a)  thia  hosmoaorphiaaa  are  unifora  with  reapect  to  r  e 
[0,T].  Thus  it  Bufficea  to  ahow 


l(A(t)  -  A(a))#lj^,  a  LI  t-a«  »/* 


for  all  #  e  H’  n  H*,  and  t,s  e  C0,T].  We  have 

0 


l(A(t)  -  A(s))Alj^,  <  t(a.(t)  -  a,(8))A,lj^,  +  i(a(t)  -  a(s))#„ij_, 

I  1 

=  (  J  '(a,(t)  -  a,(s))#,l»dx)‘/*  +  (  |  la(t)  -  a(s) l »dx] 

o  0 

I  t 

a  >/5  lAljj,  [  JlJ  axt{ff.x)d<rl*dx]  +  supla(t,x)  -  a(s,x)l  •♦'jja 
«  •  * 

t  t  t 

*  J  Jla«t(<y»’‘)'*dadx)‘'^*lt-sl*/*  +  sup  |lat(ff,x)  Ido] 


a  laig,  +  v5 


This  eatimate  implies  the  claim. 

As  an  immediate  consequence  of  Lemma  2.2  and  the  results  in  [20;  V.6 
and  V.7]  we  have: 


Proposition  2.4.  The  operators  A(t)  with  a  =  a(x,t)  e  A,  generate  an 

evolution  system  U(t,s),  0  a  s  a  t  a  t  on  L*(D).  If  f  is  Holder 

continuous  in  t  with  values  in  L’(Q),  then  the  unique  solution  u(t)  of 
(2.1),  given  by 


u(t)  =  U(t,s)A  +  J  U(t,ff)f{a)do 


satisfiss  u'c  C(tf,T;L*(Q))  and  u  c  C(8,T:H*  n  B*)  for  every  4  e  (0,T). 

0 

We  now  turn  to  the  existence  question  related  to  (P4 ).  Observe  that  A 
is  a  bounded,  cloeed  and  convex  subset  of  H'  (Q).  Therefore  it  is  also 
weakly  closed  and  weakly  sequentially  compact  in  HMQ)« 


Propoaitlon  2.5.  If  u  ^  •  or  S  ^  0  then  there  exists  a  solution  a  of 
(P,)fl»M  for  I  =  1,2. 


I 


I* 

«* 

!•* 


S«A 


•mt 


*  •*i-**4  »*>•«' 


Proof.  We  consider  the  case  i  =  2  when  a  is  a  function  of  x  and  t. 
The  cases  when  i  =  I  or  when  the  coefficient  is  a  fucntion  of  x  only 
are  treated  quite  similarly.  Let  aj  be  a  minimizing  sequence,  so  that 


J8(®j)  ■*  inf  Ji(a)  =  d 

A 


If  M  ^  *1  then  (aj  1  is  a  bounded  subset  of  .  If  Q  ^  0  then  for 
every  e  ^  0 


for  all  sufficiently  large  j.  In  either  case  (sj)  is  a  bounded  subset  of 

R*(Q).  Thus  there  exists  a  std>8equence  a  of  a.  converging  weakly  in 

Jk  ■* 

H’(Q)  to  some  a  e  H*(Q).  Since  A  is  weakly  closed  a  e  A.  Prom  (2.3) 

and  since  (a  )  is  also  a  bounded  subset  of  C(Q),  there  follows  the 
Ik 

existence  of  a  constant  k  independent  of  a  such  that 

Jk 


(2.6) 


L»(H~‘) 


By  (2.6)  there  exist  u  e  W(0,T}  and  a  subsequence  of  a  ,  again  denoted 

Jk 

by  a^  ,  such  that  u(a  )  •*  u  weakly  in  W(0,T).  Henceforth  we  suppress 
Jk  Jk 

the  index  k.  Observe  that  due  to  the  weak  convergence  of  aj  to  a  in 
H’(Q)  it  follows  that  a.  ->  a  strongly  in  C(Q).  Let  w  e  &(0,T;H*)> 

J  O 

which  is  the  space  of  all  infinitely  differentiable  functions  vanishing  at 
the  end-points  with  values  in  Let  dcoote  the  inner 

product  L'(Q)  and  take  the  limit  in 


denote  the  inner 


<Ut(aj).w>L,(Q)  =  -<aju(ej)„w.>j^,^Qj  + 

to  obtain 

(2.7)  ^'^’**^L*(Q)  ”  ^^***^L*(Q)' 

Since  l>(0,T:H*)  is  dense  in  W(0,T)  [19,  p.  11],  (2.7)  holds  for  all 

w  <  W(0»T).  Moreover  W(0,T)  embeds  continuously  into  C(0,T;L*)  and 
therefore  u(0)  =  #.  Together  with  (2.6)  and  uniqueness  of  the  solution  to 
(2.1)  in  W(0,T)  we  have  u  =  u(a).  Finally  by  the  weak  lower 


Mb;: 

i5: 

m 

7* 
% 

I 

I 


•emicontinuity  of  the  norm 


d  =  lig  Ja(aj)  =  lig  (lu(-,T;aj)  -  **'^*(0) 


H»(Q) 


‘  iiS  .T:aj)  -  Za'La(n)  ^  ®  ^ 


‘  T:a)  -  .  Q'»'Sa(Q)). 


Therefore  a  ia  a  solution  of  (Pj)  and  the  claim  ia  verified^ 


n.  fnr  ■pAtially  dependent  diffuMon  coefficient*. 

In  this  section  wc  examine  the  case  in  which  a  ia  a  function  of  x 
only.  Thus  A  =  (a  e  H*(0):  a(x)  ^  i*,  laljj!  e  ax*d  the  regularization 
ters  in  J,  has  the  fors  Qiai*, For  a,  a  (local  or  global)  solution 
of  (Pt)^i  we  define  the  open  set 

U  =  U(a<)  =  (x  €  Q:  a,(x)  ^  v) , 

and  assume  throughout  that  U  is  not  empty.  Recall  that  U(a,)  can  he 
written  as  a  countable  union  of  open  intervals.  Throughout  the  first  part  of 
this  section  we  consider  the  case  of  the  regularized  fit-to-data  criterion, 
i.e.  fl  ^  0.  We  have  the  following  results: 

Theorem  3.1.  Let  ai  bo  a  solution  of  (Pi)®  with  fl  ^  0  (and  v  ^  ^  • 
arbitrary).  Then  tit  e  H>(U).  If  aoreover  ♦  e  H*  and  f  e  L*(L*),  then 
a,  e  H»(U). 


Theorem  3.2.  Let  a^  bo  a  solution  of  (Pa)®  with  0  ^  0  (^  arbitrary). 
Then  a,  e  W**‘(U).  If  ♦  c  H*  and  f  e  L*(L*)  then  a,  t  H»(U).  If 
aoreover  e  H*  then  aj  e  H’(U). 

In  addition  to  these  local  regularity  properties  of  a^  we  obtain  the 
following  global  regularity  property  for  at*. 

Theorew  3.3.  Let  fl  ^  0  and  let  ai  be  a  solution  of  (Pf)®,  i  =  1,2.  Then 
a^  is  of  bounded  variation  and  in  particular  a  e  H''*'*(0,1)  for  every  s  e 
C0,i/*). 

The  case  0=0  is  considered  at  the  end  of  this  section. 

To  prove  these  theorems  we  prepare  several  preliminaries.  First  the 
Pr4chet  derivative  of  the  parameter-to-output  mappings  4t:  ^  H* 

W(0,T),  given  by  ft  (a)  =  u(a)  and  0  «=  H*  L*(Q)  given  by  ♦i(a)  = 
u(*,T*,a)  is  characterized.  For  a  c  A  and  h  <  H*  consider: 


Vt  =  (av,),  +  (hu,(a)),  in  Q, 


(3.1)  v(x,0)  =  0  in  n 

v(0,t)  =  v(l,t)  =  0. 

Since  (hu,(a)),  e  L*(0,T;H~‘)  Proposition  2.1  implies  the  existence  of  a 

unique  solution  v  e  W(0,T)  of  (3.1).  If  aoreover  ♦  e  H‘  and  f  c  L*(L*) 

0 

then  by  Proposition  2.3  u<a)  e  L*(H*)  and  one  can  conclude  that 
{hu(a)*)*  e  L*(L*).  Applying  Proposition  2.3  once  again  wo  find  that  v  e 
L*(H*)  with  €  L*(L*).  The  Frechet  derivative  of  ♦)  in  a  in 

direction  h  is  denoted  by  it*(a;h). 


Proposition  3.1.  a)  The  Frechet  derivative  of  ♦,  in  a  e  A  in  direction 
h  €  H‘  is  given  by  the  unique  solution  v  of  (3.1),  i.e.  4i(a;h)  =  v. 
b)  Similarly  W(a;h)  =  v(-,T). 


Proof.  We  shall  use  the  implicit  function  theorem  [6,  p.  115]  and  define 
F:  H‘(Q)  *  W(0,T)  L*(0,T;H-‘)  >  L*(Q)  by 

F(a,u)  =  (ut  -  (au,),  -  f,u(0)  -  ♦). 

One  can  easily  show  that  F  is  continuous  and  that  Fu(a,u),  a  e  A,  is 
surjective  from  W(0,T)  to  L*(0,T;H"‘)  *  L*.  Moreover  the  partial  Frechet 
derivative  of  F  with  respect  to  a  in  direction  h  is  given  by 

F,(a,u)(h)  =  {(hu,)„,0).  Clearly  (a,u)  -»  F,(a,u)  from  H‘ (0)  »  W(0,T)  to 

£(H*,L*(H~')  X  L’(0))  is  continuous.  Here  £(X,Y)  denotes  the  set  of  all 
bounded  linear  operators  from  X  to  Y.  Thus  by  the  implicit  function 
theorem  a  4i(a)  =  u(a)  is  Frechet  differentiable  from  H*  to  W(0,T). 
Its  value  can  now  be  calculated  directly  from  (2.1)  and  is  found  to  be 

characterised  by  (3.1).  As  for  the  Frechet  derivative  of  #3  observe  that 

u  ■*  u(*,T)  from  W(0,T)  to  L*(Q)  is  a  bounded  linear  mapping.  Therefore 
W(a!h)  =  v(*,T).  From  Proposition  3.1  it  follows  that  JJa),  i  =  1,2  are 
Frechet  differentiable  with  respect  to  a.  The  Frichet  derivative  of  J|  at 
a  c  A  in  direction  h  e  H*  is  given  by 

Jl(s;h)  =  2<u(a)  -  ♦  2fl<a,h>g, 


and 


Ji(a;h)  =  2<u(a)(-,T)  -  *a» Va>j;^a^Qj  +  2B<a,h>„, 


H‘(Q)- 


Integration  by  parta  and  uae  of  (3.1)  yields  an  alternative  representation 
for  j/(a;h). 


Proposition  3.2.  Let  a  e  A  and  h  e  HMO).  Then 


r 

Jl{a:h)  =  -2  J  J  u,(a)p,(a)dt  h  dx  +  2fl<a,h>gi , 


0  o 


and 


Ji(a;h)  =  -2  I  J  u,(a)q,(a)dt  h  dx  +  2fl<a,h>g, , 
Q  « 


I 

Lmb8_3j;J^  Let  w(x)  =  J  Ux(x,t)rx(x, t)dt. 

o 

a)  If  u  e  W(0,T)  and  r  «  W(0,T)  then  w  e  L‘. 

b)  If  u  €  W(0,T)  and  r  e  L»(H»)  then  w  e  L*. 

c)  If  u,  €  C(L*)  and  r  e  L*(C*)  then  w  e  L*. 

d)  If  u  €  L*(H*)  and  r  €  L*(H*)  then  w  e  H‘. 

3.2.  Throughout  this  paper  and  specifically  in  the  proof  of  Lemma 
3.1  we  use  the  following  fact  which  is  a  special  case  of  [8,  Theorem  17;  p. 
198].  If  P  e  L'(0,T;LP(Q))i  1  ^  p  e  •,  then  there  exists  an  integrable  (with 
respect  to  the  product  Lebesgue  measure  on  [0,T]  x  Q)  function  f, 
uniquely  determined  up  to  measure  zero,  such  that  f(t, ‘)  =  P(t)  for  almost 
every  t  e  [0,T].  Moreover  f(',x)  is  integrable  for  a.e.  x  e  Q  and 

1  I 

J  f(t,x)dx  as  a  function  of  x  equals  the  element  J  F(t)dt  of  LP(Q). 


y* 

V 


Proof  of  Lemma  3.1.  Part  a)  is  obvious.  To  verify  b)  observe  that 


I  I  *1  I 

JlJ  u,r*dtl*dx  <  J[J(u*l*dt  •  J<r, l*dtjdx 


*  s^  Jlr.l«dt  4  const  '"'wCO.T) ' 

xeu  0 


To  verify  c)  we  use  Minkowski's  integral  inequality  [21,  p.  271]  in  the  first 


step: 


II  II 

[JlJ  u„r*dtl  *dx]  J[Jlu,r,l  *dx|  *'^*dt 

o  d  o  o 

T 

-  Jlr*(-,t)l^  lu,(-.t)lj^,dt  -  '«-'li(c«)'“*'c(L»)' 


Part  d)  is  left  to  the  reader. 


litMl  ?t2i  Let  I  be  an  open  interval  in  Q,  d  ^  0  and  a  c  HMD)*  If 
further  f  c  L‘  (or  L*  or  HM  and 


IWH  WllW  WWlWl.W.lli  W  WmiR.'  WUUiMJUMUUUWWm  V  V  M,<m Ml  Ml  PUI WVH  Wlfl  UTVT  VW).' 


I  [d(a,hx  *  ah)  +  fbldx  =  0  for  all  h  € 
i 

then  -a,*  +  a  =  -  ^  •^, 
and 

a  e  W’.'c'i)  (or  or  H*(I)). 


Proof.  By  aaaunption 

X  X 

J  (aa,  -6|  ads-J  f)h,dx  =  0 

I  0  0 

X 

for  all  h  €  H‘(I).  Therefore  a,  =  J  (a  +  j)*!®  const,  and  therefore 
®x»  ~  ^  renainlng  assertion  now  follo%«s  easily. 

The  proof  of  Theorems  3.1  -  3.3  employs  a  Lagrange  multiplier 
formulation  of  (P^).  As  another  preparatory  result  we  therefore  discuss  the 
regularity  condition  for  the  constraint  set  A.  Let  Gi.’  H^(Q)  •*  R  and  Gj; 
H‘{D)  "*  C(0)  be  defined  by 

G|(a)  =  iai*,  -  M*  and  G,{a)  =  v  -  a, 

and  put  G(a)  =  (G| (a),G](a))  where  G:  H‘  -»  R  x  C.  We  have 

G|(a:h)  =  2<a,h>g,  and  Gj(a;h)  =  -h. 

By  we  denote  the  cone  of  nonnegative  functions  in  C|  by  Rf  the 

nonnegative  reala  and  C%  stands  for  the  dual  cone  of  C4.,  i.e.: 

C|  =  (X*  e  C*;  X*(y)  a  0  for  all  f  e  C+). 

Observe  that  A  =  (a  e  H‘:  -G(a)  e  R+  x  C+>  =  (a  c  H*:  G(a)  ®  0). 


Every  admissible  point  a  c  A  satisfies  the  regular  point 
condition.  That  is 


(3.4)  0  €  inf  {G(a)  +  range  G*(a;H*)  H*.  *  C+). 


Proof.  This  reault  ia  verified  in  [IS].  For  the  aake  of  completeneaa  and 
since  the  proof  is  short  we  repeat  it  hero.  We  will  show  that  every  (r.V') 
is  contained  in  the  set  appearing  in  (3.4)  provided  that  Irl^  +  Ifl^  ^  t  = 


"^“(jTal — ’5  Note  that  (3.4) 


4lalj^,’2 


is  equivalent  to 


(3.5)  0  €  infdaljj,  -  +  2<a,h>jji  +  R4.,  i/  -  a  -  h  +  C+:  he  H->. 


Lei  f  -  f  ~  min  f  and  note  that  V  e  C4..  We  put  h  =  t/  -  a  -  min  f  e  H‘. 
Then  y'  =  i/-a-hf9,  which  ia  of  form  of  the  elements  in  the  second 
comijonent  of  (3.5).  Next  we  need  to  find  p  e  tt+  such  that 


(3.6)  r  =  lal',  -  p*  +  2<a,h>jj,  +  p. 

This  is  cK}uivalent  to  r  =  -p*  -  '«'ga  2(a,v)j^,  -  2(a,ain  f)j^,  +  p.  But 

2(a,i/)j^,  -  2(a.Bin  f)j^,-lal  *i-m*  <  lal  *a+v*-lalgi-p*-2(a,min  f)j^a 
*  v’  -  p*  +  2lal^,e  ^  ^(v*  -  p*). 

Therefore 


p  *  r  +  |(p*  -  V*)  »  |(p*  -  V*)  -  Irl  *  ^(p*  -  V*)  -  e  k  0. 
Thus  p  in  (3.4)  ia  nonnegative  and  the  lemma  is  verified. 


Proof  of  Thaoram  3.1.  We  only  consider  the  case  p  ^  The  case  p  =  • 
requires  no  esaential  changes  from  the  preaant  proof. 


By  the  Lagrange  multiplier  theorem  (see  e.g.  [24])  there  ezista  a  pair 
(XiaAf)  c  14  X  C|  such  that 

(3.7)  j;(a;h)  +  A|G;(a:h)  ♦  <A!,G;(a;h)>^,  +  B<a,h>g,  =  0, 

for  all  h  c  H*, 


IIWWVU  tfWtnWWlWIlWJWV.tf^'V'JWWinniWH^VVTrLniWWW^^  yjir,- 


(3.8)  Ai(l®l*,  -  /la)  =  0  and  <xf,i/  -  a>^  =  0. 

Here  denotes  the  duality  pairing  on  C  and  we  dropped  the  index 

1 

in  the  notation  for  ai .  Recall  that  Xf  can  be  expressed  as  Xf  =  J  dXa» 

O 

where  Xa  is  a  nondecreasing  function  of  bounded  variation.  From  the 
mean  value  theorem  for  Stieltjea  integrals  and  (3.8)  it  follows  that  Xa  is 
constant  on  open  subintervala  of  U  =  U(a).  Using  Propoaition  3.2  it  can  be 
seen  that  (3.7)  ia  equivalent  to 

(3.9)  -2(w,h)j^a  +  (2Xi+fl)<a,h>jja  +  <xf,h>^  =  0 

for  all  h  e  H*.  where  w(x)  =  (J  Ux(a)Px(a)dt) (x) .  Let  I  be  an  open 

o 

interval  in  U  and  h  e  (h  e  H*:  h  =  0  on  Q  XI).  Then 

(3.10)  ~2(W|  h)  j^a  ^  J  J  (  2Xi+B)  <a(  h^jji  ^  J  J  -  0. 

By  Remark  3.1  a)  and  Lemma  3.1  b)  we  have  w  e  L‘(0)  and  under  the 
assumption  #  e  and  f  e  L’(L’)  it  follows  from  Remark  3.1  a)  and 
Lemma  3.1  d)  that  w  e  H*.  Since  0^0  by  assumption,  we  can  use  Lemma 
3.2  with  A  =  fl  -f  2X,  ^0  to  obtain 

(3.11)  -a,,  +  a  =  w  on  I 

and  a  e  H*(I)  respectively  a  c 

Since  U  is  an  open  subset  of  Q  in  dimension  one,  it  can  be 
expressed  as  a  countable  union  of  disjoint  open  intervals.  On  each  of  these 
intervals  (3.11)  holds  and  therefore  a  c  H*(U),  respectively  a  c  H’(U)  if 
the  stronger  regularity  assumptione  hold  for  f  and  #. 


Proof  of  Theorem  3.2.  The  proof  is  quite  analogous  to  the  one  of  Theorem 

T 

3.1.  Here  we  take  w(x)  =  (J  Ua(aa)qK(aa)dt) (x)  and  by  Ramark  3.1  b  and 

0 

Lemma  3.1  a)  we  find  w  c  L*.  Again  we  have  (3.11)  and  therefore  a  c 


and  consequently  a  e  W*i*(U).  If  in  addition  ♦  e  H*  and  f  e 

0 

L*(L*)i  then  u,  e  C(L*)  and  q  «  LMC‘)  and  thus  w  e  L*  by  Lemma 
3.1  c).  (fader  the  additional  assueption  e  H‘  we  find  w  e  H*  by 

Remerk  3.1  c  and  Lemma  3.1  d).  Using  these  facts  in  (3.11)  we  arrive  at  the 
desired  result. 


Proof  of  Theorem  3.3.  Consider  (3.9)  with  w  defined  as  in  the  proof  of 
Theorem  3.1  if  i  =  1  or  as  in  the  proof  of  Theorem  3.2  if  i  z  2.  In  either 
case  w  c  LMQ).  Thus  from  (3.9)  the  mapping  a,  ((ai)x,h,)  can  be 

1 

extended  to  a  bounded  linear  functional  on  C,  i.e.  (ax,hx)  =  |  h  for 

0 

all  h  e  and  some  function  i>  of  bounded  variation.  For  all  h  e 

o 

we  thus  have 

1 

(a,,h,)  =  -  J  f  h,dx 
0 

which  implies  a*  =  -f  +  C  for  a  constant  C.  Therefore  a  is  of  bounded 
variation  and  the  proof  is  finished. 

Remark  3.2.  It  is  interesting  to  observe  that  in  the  case  A|  =8  =  0  the 
regularity  of  Aj  (which  is  used  for  the  ■daaekeee  integral  representation 
of  a})  is  determined  by  w.  In  fact,  in  this  case  w  =  Aj  and  Aa  e  W*  >  ‘ 
by  (3.9)  for  i  =  1  or  2. 

Now  let  us  consider  the  case  8  =  0.  From  (3.9)  -  (3.11)  it  is  obvious 

that  the  aame  regularity  properties  can  be  obtained  for  the  solutions  ai 

of  the  unregulariied  problems  (P^)*,  i  =  1,2,  if  only  A|  «  0.  By  (3.8) 
this  la  turm  implies  that  the  norm  constraint  is  active,  i.e. 

We  recall  that  U(a<)  is  assumed  to  be  nonempty. 

PropoelUon  3.3,  Let  8  =  0,  m  ^  •  and  let  a,  be  a  solution  of  (P,)*.  If 

for  some  open  interval  1  ^  q 

T 

(3.12)  J  Ux(&i)pB(a,)dt  «  0  in  L*(I)  respectively 


I 

J  u«(®j)q«(«a)dt  *  0  in  L*(I), 


then  the  Lagrange  multiplier  associated  with  the  norm  constraint  satisfies 
A|  ^  0  and  the  regularity  results  of  Theorem  3.1  -  3.3  hold.  Moreover 
=  M  in  this  case. 


Proof.  Assume  that  Xi  =0  and  choose  I  as  in  the  statement  of  the 
proposition.  Then  by  (3.10)  we  have  w  =  0  on  I  which  contradicts  (3.12) 
and  ends  the  proof. 


It  appears  to  be  infeasible  to  give  conditions  in  terms  of  (2.1)  and  the 
equations  for  p  or  q  which  guarantee  (3.12).  The  following  proposition 
gives  a  necessary  condition  for  (3.12)  to  hold  and  thus  for  X|  ^  0. 


Proposition  3.4.  Let  0:0  and  u  ^  It  X|  ^  0  then  u(a4 )  cannot 
coincide  with  Zf  as  a  function  in  L*(L*)  respectively  L’(Q). 


Proof.  Let  I  c  u  be  an  open  interval  and  assume  that  uia^)  =  Z|,  for 
i  =  1  or  2.  Then  p  =  0  respectively  q  =  0.  Therefore  (3.12)  cannot 


hold.  By  (3.10)  =  0  for  all  h  e  R*(I)  and  thus  X|  =  0. 


ss 


4.  Re. 


I 

1 

I 


f  for  aT<«tt»lly  wnd  teaporally  dependent  diffuaion  coefficienta. 

In  thia  aaction  we  consider  the  caae  in  which  the  coefficient  a  is  a 
function  of  both  x  and  t.  Hence  we  take  the  admissible  set  of  parameters 
to  be 


A  =  (a  €  H*(Q);  a(x,t)  *  i'  ^  0,  'aljj,  s 

and  the  regularization  tern  in  J(  is  of  the  fora  ^*^*H*(q>*  ^ 

(local  or  global)  solution  of  (P|)^  we  define  the  open  set 

U  =  U(a<)  =  ((x,t)  €  fl:  ai(x,t)  ^  v) . 


Throughout  thia  section  it  ia  assume  that  U  is  not  empty.  As  in  the  case 
of  section  3,  the  optimal  parameter  may  satisfy  special  regularity  properties 
either  due  to  the  regularization  term  or  due  to  the  norm  bound.  Here  we 
state  both  cases  simultaneously.  We  shall  write  0|  q,  if  and  Oa 

are  open  subsets  of  a  Euclidean  space  and  if  fl|  is  compact  with  Qi  c  Oa*  < 

I 

The  Lagrange  multiplier  associated  with  the  norm  bound  is  again  denoted 
by  X,. 


Theorem  4. 1.  Assume  that  li  ^  •  and  X|  *  0  or  that  B  ^  0  and  let  at 

be  a  solution  of  (P,)^.  If  U  ««  U,  f  c  L»{Q)  and  ♦  e  H‘ ,  then  a  e  H*(U). 

0 


Theorem  4.2.  Assume  that  n  ^  •  and  X|  «  0  or  that  0^0  and  let  ha 
be  a  solution  of  (Pj)^.  If  U  cc  u,  f  «  L’(Q),  #  e  and  Zj  e  then 

0  O 

a  c  H«(U). 


Theorem  4.3.  Assume  that  u  ^  *  mnd  Xt  «  0  or  that  0^0  and  let  a^ 
be  a  solution  of  (P|)^>  i  =  1  or  2.  If  Q  cc  q  and  Q  satisfies  the 
uniform  cone  property,  then  Si  e  H’'^*(Q)  for  any  s  c  (0,1). 


In  section  5  we  shall  specify  conditions  which  guarantee  that  Xi  *  0. 


The  verification  of  the  above  theoremsi  proceeds  along  the  same  lines 
as  that  presented  in  section  3.  We  indicate  the  salient  features.  First  the 


Pr^het  derivatives  of  J),  1  =  1  or  2,  are  specified. 


Proposition  4.1.  Lot  a  e  A  and  h  £  H*{Q).  Then 

T 

Jl(a:h)  =  -2  j  I  u,(a)p*(a)h  dtdx  +  2B<a,h>jja^Qj , 
n  o 


and 


Ji(a:h)  =  -2  I  I  u,(a)q,(a)h  dtdx  +  2fl<a,h>jja^Qj , 

0  o 

where  p  e  W(0,T)  and  q  e  W<0,T)  are  given  in  (3.2)  and  (3.3). 

Remark  4.1.  Since  f  £  L*(H”*)  and  #  e  L*  it  follows  that  u  e  W(0,T) 
and  that  u,  e  L^(L^).  By  Proposition  2.3  this  further  implies  that  p  e 
L»(H*)  and  pt  e  L»(L*). 


Lemma  4.1.  a)  Under  the  standard  assumption  that  f  e  L’(H~')>  #  e  L*,  z, 
£  L*(Q)  and  Zj  e  L’fQ)  it  follows  that  u^p,  e  LMQ)  and  u,p,  e  L'(Q) 
and  thus  both  these  functions  can  be  identified  with  elements  in  the  dual 
space  (H*'^V)*  of  H‘'*’y(Q)  for  every  y  ^  0. 

b)  If  f  e  L*(L*)  and  #  c  then  e  L’(L’). 

c)  If  f  c  L*(L‘),  #  e  and  Zj  e  then  u^q,  e  L*(L*). 


Proof.  Obviously  u^Px  and  u^q.  belong  to  L'  (Q).  Moreover,  if  f  e 
L*(Q)  then,  since  H*'^'>'(Q)  embeds  continuously  in  L*(Q)  for  y  ^  0,  we 

see  that 

=^^JV?y  '♦'g.+y 

for  aome  constant  c  independent  of  ft  L^(Q).  Elere  denotes 

the  duality  pairing  between  and  with  L*  as  a  pivot 


■pace*  To  verify  b)  we  observe  that  by  Proposition  2.3  we  find  that  u  e 
C(H‘)»  Moreover  p  e  L*(H*)  by  Remark  4.1  and  thus  it  easily  follows  that 
u,Pa  e  L*(L*).  Finally  c)  is  iaplied  by  the  fact  that  u  c  C(H‘)  and  q  c 

To  determine  the  regularity  properties  of  the  solutions  of  (P^)^ 

we  again  use  a  Lagrange  formulation.  First  the  constraints  that  are 
involved  in  defining  the  set  of  admissible  parameters  are  characterized  by 
a  mapping  G  =  (Gi,Ga):  H*  -*  R  x  given  by 

Gi(a)  =  'a'ga(Q)  -  M*.  Ga(a)  =  1/  -  a. 

Let  =  (a  €  a(x,t)  »  0  on  Q>.  For  a  e  we  write  a  <  0 

if  -a  «  Observe  that  a  e  A  if  and  only  if  G(n)  a  0.  Since 

y  a.  0,  embeds  continuously  into  C|  one  can  show  writh  an  argument  similar 
to  the  one  of  Lemma  3.3  that  every  point  a  e  A  satisfiea  the  regular  point 
condition 


0  e  int(G(a)  +  range  G’{a;H*)  +  Rf  x 

Thus  for  the  problems  (Pj)®,  i  =  1  or  Z,  there  exists  a  Lagrange 
aultiplier  pair  (X,,a|)  e  R+  x  satisfying 

T 

(4.1)  -2  J  J  u,(ai)p,(a,)h  dxdt  +  2<Xt+Q><a, ,h>g,^Qj  +  <x!,h>j^y  =  0, 

o  Q 


respectively 

T 

(4.2)  -2  J  J  u,(a,)q,(a,)h  dxdt  +  2<X,+B><a,,h>g,^Qj  +  <xf,h>j^y  =  0 

0  Q 

for  all  he  H*(Q).  Moreover  the  complementary  condition 

(4.3)  A»(»a»lJ,  -  /i")  =  0,  <A,,y  -  a,>^^y  =  0, 


respectively 


(4.4)  Xx(la,l',  -  M*)  =  0.  <X,.y  -  =  0. 

holds.  Here  =  (g  «  »  0  for  all  a  e  H‘+V) .  Next 

we  analyse  furhter  the  variational  equations  (4,1)  and  (4.2).  Henceforth  the 
assumption  X|  Q  ^  0  of  Theorems  4.1  -  4.3  will  be  used.  The  bilinar 
form 

B{a,h)  =  2(X,+fl)<a,h>jj, 
is  H’-coercive  and  H’-bounded,  i.e. 


B(a,a)  a  2(X|-t-Q)  lal|]2 


B(a,h)  -  2(X,+a)iaijj,ihlg,, 

Thus  by  the  4^u-Milgram  theorem  there  exists  a  bounded  linear  operator 
L:  H*(Q)*  HMQ)  satisfying 

(4.5)  B(L(p),h)  =  <p,h>^jjj^S  for  all  h  e  H*(Q) 


(4.6)  '^<^)'h*(Q)  ^  2(X,+B)  '^'h>(Q)*' 

Here  ga  denotes  the  duality  pairing  between  H*(Q)  and  H*(Q)*. 

Hence,  by  restriction,  for  any  Q  cc  q  ^  find 


•I.(p)lHa(Q)  *  2(X,+0)  '^'h»(Q)*’ 


and  L  can  be  considered  as  a  bounded  linear  operator  from  H*(Q)*  to 
H*(Q).  We  shall  write 

L  c  I(H»(Q)*,H*(Q)). 

where  £(X,Y)  denotes  the  space  of  all  continuous  linear  operators  between 


the  Banach  spaces  X  and  Y.  Applying  results  on  the  interior  regularity 
of  solutions  of  variational  elliptic  boundary  value  problems  we  further 
obtain  that  for  p  e  L*(Q)  and  Q  Q  the  solution  L(p)  of  (4.5) 
satisfies  constant  C  is  independent  of 

p  e  L’{Q)  [1,  p.  51-81;  23,  p.  299-302].  We  thus  obtain 

(4.8)  L  e  I(L^Q).H*(Q)). 

We  next  utilize  results  from  interpolation  theory.  Since  H^(Q)  is 
continuously  and  densely  imbedded  in  L*(Q)  it  follows  that  (taking  L^(Q) 
as  a  pivot  apace)  L^(Q)  can  be  identified  with  a  subset  of  (H^)*. 
Subsequently  from  (4.7),  (4.8)  and  the  interpolation  theorem  we  find,  using 
the  notation  from  [19], 

(4.9)  L  €  I([L*(Q).H*(Q)*]g,(H-(Q).H»(5)]g) 

for  every  S  e  (0,1).  For  interpolation  between  dual  spaces  it  is  known  [19, 
p.  29]  that 

[L*(Q).H»(Q)*]g  =  ([HM0).L»(Q)],-9)*. 
with  equivalent  norma,  and  further  [19,  p.  40] 

([H*(Q),L>«l)]»-g)*  =  H*«(Q)*. 

This  implies  [L*  (Q),H*  (Q)*  ]g  =  H*®(Q)*.  If  Q  satisfies  the  uniform  cone 
property,  then  the  Calderon  Zygmund  extension  theorem  [23,  p.  100]  is 
applicable  to  the  domairi  Q  and  the  results  in  [19,  p.  40-43]  can  be  used 
to  obtain 

[H*(5),H«(Q)]g  =  H»(»-®)(Q). 

Using  these  facts  in  (4.9)  we  arrive  at 

( 

i 

(4.10)  L  €  X(H»®(Q)*,e*(»-«)(Q)).  I 

I 

for  any  9  c  (0,1)). 


Proof  of  Theorem  4.1.  Let  a,  be  a  solution  of  (P,)^  and  assume  that  U 
cc  u.  Then  there  exists  a  bounded  dooiain  V  such  that  0  V  cc  \j. 


Moreover  there  exiete  5^0  such  that  ai(x,t)  ^  1/  A  on  V.  Let  h  e 
H*  with  coapact  support  in  V+  satisfy  'hlj^,  *  fl.  Then  *  h  +  i'  ■  a,  s  0 
and  by  the  complementarity  condition  (4.3) 

<A*,*h>i+y  =  <X*,*  h  +  v  -  ai>i.t-y  s  0. 

Therefore  it  follows  that 

(4.11)  <A*,h>i+y  =  0  for  all  h  c  H*  with  cosipact  support  in  V. 

We  use  this  fact  in  (4.1)  and  find 

( Ai+B)  <ai  -  0 

for  every  h  e  with  compact  support  in  V.  From  Lemma  4.1  we  have 
UxP„  e  L'(L*).  The  assumptions  in  Theorem  4.1  guarantee  that  A|  4-  0  ^  0. 
Thus  from  (4.5)  and  (4.8)  with  Q,  Q  replaced  by  V,  U  we  obtain  ai  e 
H*(U).  See  also  [23,  p.  73,  299].  This  ends  the  proof. 

Proof  of  Theorem  4.2.  We  only  replace  (4.1),  (4.3)  by  (4.2),  (4.4)  in  the 
proof  of  Theorem  4.1  to  obtain  the  desired  result. 

Proof  of  Theorem  4.3.  From  Lemma  4.1  we  observe  that  u,p,  and  u,q„ 
can  be  identified  with  elements  in  H‘*y(Q)*.  Furthermore  A*  e  H‘'*'‘>'(Q)* 
as  well.  Therefore  we  have 

Pi  :=  -Af  +  Uk(*»)Px(«i)  «  H‘'^V(Q)* 

and 

Pa  :=  -Af  +  u,(a,)q,(aa)  e  B*'^y(Q)*‘ 

By  assumption  A|  -f  Q  ^0.  Thus  from  (4.1),  reai>ectively  (4.3),  (4.5),  and 
(4.10)  it  follows  that  L(P|)  e  H’")'  for  every  y  e  (0,1).  Setting  a  =  1  >y 
this  implies  the  result. 


In  this  section  we  give  conditions,  which  guarantee  that  A,  ^  0.  In 
view  of  (4.3),  (4.4)  these  conditions  imply  that  the  norm  constraint  is  active 
and  they  constitute  a  specific  case,  where  the  strict  complementarity 
conditions  holds: 

“  y*)  =  0,  i  =  1  or  2,  without  both  factors  being  zero 
simultaneously.  We  shall  see  that  the  problem  of  nontriviality  of  A,  is 
related  to  the  question  of  attainability  of  z«.  Here  Z(  is  called  attainable, 
if  Zi  €  Vi  z  (z:  z  z  u(a)  in  W(0,T),  a  e  A),  respectively  Zj  e  Vj  =  (z:  z 
=  u(’,T;a)  in  L’(a),  a  e  A),  where  u(a)  is  the  solution  of  (2.1). 

We  point  out  some  of  the  consequences  of  the  nontriviality  of  A,. 
First,  we  have  already  seen  in  sections  3  and  4  that  At  ^  0  implies 
regularity  properties  of  the  solutions  a^  of  (P« ).  Moreover,  for  elliptic 
estimation  problems  At  ^0  guarantees  stability  of  the  solution  of  the 
output-least-squares  formulation  of  the  estimation  problem  with  respect  to 
perturbation  of  z  or  A  C7]  and  it  guarantees  augmentability  as  used, 
for  example,  in  the  augmented  Lagrangian  approach  to  estimation  problems 
that  we  described  in  CIO]. 

Recall  from  (4.1)  that  for  0  =  0 

(5.1)  — 2^u„p„ ,  hA|^2  +  2A I  <at ,  hA||]  +  ^Af.hAj^^  =  0 

for  all  h  e  and  that 

U  =  U(ai)  =  ((x,t)  €  Q;  ajCx.t)  ^  v) . 

For  the  unregularized  fit-to-data  criterion  J|  we  have  the  following 
result  which  guarantees  A,  ^  0.  By  a  rectangle  V  in  Q  we  mean  a  set 
of  the  form  [xt,xj]  x  [t|,ta]  <=  Q  with  0  <  X|  ^  Za  «  1  and  0  e  ti  ^  t, 
<  T. 


Theorem  5.1.  Let  B  =  0,  f  e  L’(L*),  #  e  H‘,  Zj  be  Holder  cootinuous  with 
values  in  L*(Q)  and  assume  that  ai  is  any  (local  or  global)  solution  of 
(Pi)  satisfjring  ki  0  v  and 

(a)  for  every  rectangle  V  c  u,  f  is  not  (a.e.)  only  a  function  of  t  on 
V, 


(b)  there  exists  at  least  one  rectangle  V  c  u  such  that  u(at)  -  Z|  is 

not  (a.e.)  only  a  function  of  t  on  V.  Then  Xi  3^  0. 

Before  we  give  the  proof  let  us  further  interpret  the  result.  First 

observe  that  the  conclusion  Xt  ^0  depends  on  u  only  in  as  far  as  ai 

may  depend  on  u.  The  solution  kt  of  (P,)  enters  into  the  assumptions 

only  in  (b)  and  through  the  requirement  that  ki  P  n  (which  is  equivalent 

to  U(ai)  *  0).  As  for  (b)  we  point  out  that  due  to  the  assumptions  on  f 

and  k  every  solution  u(a),  a  e  A,  of  (2.1)  is  continuous  on  Q  with 

u(0,t;a)  =  u(l,t;a)  =  0.  In  the  specific  case  that  z,  is  constructed  from 

piecewise  constant  interpolation  of  point  data,  (b)  holds  for  every  a  e  A 

unless  Zt  ■  0.  We  observe  also  that  under  the  smoothness  assumptions 

and  assumption  (a)  of  Theorem  5.1,  (P,)  with  the  norm  constraint  eliminated, 

cannot  have  a  solution  at,  which  is  not  identically  v  and  which  satisfies 

(b).  In  fact,  if  this  were  the  case  at  were  also  a  solution  of  the  norm- 

constraint  problem  with  p  =  2tatlg2  and  Theorem  5.1  were  applicable  and 

gave  Xt  V  0.  By  the  complementarity  condition  (4.3)  this  would  imply  p  = 

«» 

”bich  is,  of  course,  impossible. 


Proof  of  Theorem  5.1.  Let  us  assume  that  X|  ^  0.  Then  by  (5.1)  we  have 

(5.2)  2<u,p,,h>j^,^Qj  =  <A!,h>j^y  for  all  h  e  H». 

Thus  h  ^  extended  to  a  continuous  linear  functional  F 

on  L*(0)  with  F(h)  =  <2uxP*,h>j^,^Qj .  Here  we  used  the  fact  that  u  e 
C(H*)  and  p  €  L*(H*  n  H*)  by  Proposition  2.3,  so  that  u,p,  e  L*(L*). 
Since  xf  is  in  the  positive  dual  cone  we  find  that 

*  0 

for  all  h  c  n  c*.  But  o  C*  is  dense  in  L*  (this  can  be  seen 

with  a  mollifier  argument  [ADAMS,  p.  30])  and  thus  (u,p,,h)^,  »  0  for  all 
h  <  L*.  This  implies 

(5.3)  UkPk  ^0  a.e.  on  Q. 

Moreover,  since  <xf,v  -  s>,+y  =  2(u*p*,v  - 


=  0  we  find  that 


Next  observe  that  t  ■*  u(t,‘)  from  [0,T]  to  L*(Q)  is  Holder 

continuous  with  exponent  1/2.  In  fact, 

t  T 

lu(t)  -  u(s)i  *  J  iut(<y)ldff  <  it-el‘/«  (J  iutl*da)*/», 

•  0 

where  the  exisitence  of  the  last  integral  follows  from  Proposition  2.3.  Since 
Z|  is  Hblder  continuous  with  values  in  L*(Q),  it  follows  that  u  -  Zi  is 
Holder  continuous  with  values  in  L*(Q)  and  thus  p  e  C(0,T-tf:H*  n  h‘)  for 
every  A  e  (0,T)  by  Proposition  2.4.  Specifically  p^  is  continuous  on 

[0,T-63  a  tO,l]  for  every  6  e  (0,T). 

Now  let  us  return  to  (5.4)  and  assume  that 

(5.5)  PK(x,t)  *  0  at  some  (x,t)  e  U  «  ([0,1]  a  [0,T-4l). 

(Of  course,  it  is  assumed  that  A  is  chosen  sufficiently  small  so  that  U  o 
([0,1]  a  [0,T-A])  «  0.)  Then  there  exiats  a  rectangle  V  r  V(x,t)  ~  [xi,Xa]  x 
[ti,t]]  such  that  Px(x,t)  «  0  for  all  (x,t)  e  V.  By  (5.4)  it  follows  that 
u„(x,t)  X  0  on  V  and  thus  u,  is  a  function  of  t  only  on  V;  i.e. 
u(x,t)  X  C|(t),  where  W‘**(t,,ta;R).  Using  this  in  (2.1)  we  arrive  at 

cl(t)  =  f(t,x)  on  V. 

This  contradicts  assumption  (a)  of  the  theorem  and  therefore  (5.5)  cannot 
hold  at  any  (x,t)  e  U  o  ([0,1]  a  [0,T-A]),  where  A  e  (0,T).  We  conclude 

that  Px(x,t)  X  0  a.e.  on  U.  Let  now  V  x  [xi,xa]  a  [tt,ta]  U  be  chosen 

as  given  by  assumption  (b).  Then  p(x,t)  x  ca(t)  on  V  for  some  Ca  e 
W'**(it,ta;R).  Using  this  in  the  equation  for  p  we  arrive  at 

ci(t)  =  u(t,x)  -  z,(x,t)  on  V, 
which  contradicts  (b). 

Thus  we  have  to  revoke  the  assumption  A|  x  0  and  the  theorem  is 
proved. 


ftfirh  ffilr  observe  that  in  the  proof  of  the  previous  theorem  the  role 
of  u  suid  p  can  be  reserved,  provided  u  is  sufficiently  regular  (i.e. 
u.  is  continuous  on  [0,1]  a  [A,T],  for  any  A  c  (0,T)).  Specifically  one 
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can  show:  if  the  assumption  of  Holder-continuity  of  Z|  is  replaced  by 
Hdlder-continuity  with  values  in  L*(Q)  of  f  and  if  (a),  (b)  in  Theorem 
5.1  are  replaced  by 

(a')  there  exists  some  rectangle  V  c  u  f  jg  (a.e.)  only  a 

function  of  t  on  V, 

(b’)  for  every  recUngle  V  e  U,  u(a)  -  z,  is  not  (a.e.)  only  a  function  of 
t  on  V, 
then  Xt  ^  0. 

We  now  consider  the  case  of  the  fit-to-data  criterion  Jj.  We  have 
-2<u,q,,h>j^j^Q^  +  2X,<a,h>ga  +  <xf,h>j^y  =  0,  for  all  h  e  H». 

Theorem  5.2.  Assume  that  a*(x,t)  ^  i/  on  [0,1]  «  [0,Tl  and  that 

(a)  for  every  rectangle  V  <=  [0,1]  *  [0,T],  f  is  not  (a.e.)  a  function  of  t 
on  V  only, 

(b)  u(x,T;a)  is  not  a.e.  equal  to  Za(x)  on  Q. 

Then  Xt  ^  0. 


Proof.  Since  a*(x,t)  ^  v,  the  complementarity  condition  implies  that 

x|  z  0.  Therefore 

-  <u,q«,h>j^a^Qj  +  Xi<a,h>ga  =  0  for  *11  b  c  H*. 

Now  assume  that  X|  =  0.  Then 

(5.6)  u^q,  =  0  a.e.  oo  [0,1]  »  [0,T]. 

From  Proposition  2.4  we  find  that  q  e  C(0,T-A;H*  n  h^)  for  every  6  e 
(0,T).  Specifically  q„  is  continuous  on  [0,1]  x  [0,T).  If  qx(x,  t)  »  0 
for  some  (x,l)  e  [0,1]  ■  [0,T)  then  u,  =  0  a.e.  on  some  rectangle  V  ^ 
[0,1]  X  [0,T].  This  contradicts  (a)  and  thus  q,(x,t)  =  0  on  [0,1]  x  [0,T). 
Now  wo  uos  the  equation  for  q  and  find  that 

IT  T  1 

~  f  f  *****  'll  ^**1*)k^  *****  =  0 

•  •  0  0 
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